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ABSTRACT 

In this thesis the formation and properties of a polymer gel on and at a surface 
are investigated. The gel under investigation is defined as a three-dimensional net- 
work of macromolecules that form permanent links with one another and also with 
confining planar surfaces. The precise location of the crosslinks on the wall or on 
another macromolecule is not known prior to linking, and will differ from sample to 
sample. However, once the crosslinks are formed, they are assumed to be perma- 
nent. This random linking is the source of the disorder in the system, over which 
a quenched average has to be taken. An existing model [H] of network formation, 
with polymer-polymer crosslinks, is extended to incorporate a surface and polymer- 
surface crosslinks. Within the framework of replica theory, statistical averages and 
physical properties of the system are calculated by means of a variational approach. 
Macroscopic information, in terms of the free energy of deformation, is obtained by 
using two different potentials to simulate the crosslinks mathematically. 
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Chapter 1 

Introduction 



Polymers in restricted geometries exhibit very specific and interesting properties. The widespread 
interest in polymers in confined media, is proof of the importance of understanding these prop- 
erties, which play a vital role in many industrial applications, but also in nature 1 . Fundamental 
theoretical physics research provides the basis for investigating and predicting properties of these 
systems. Problems related to this field are diverse, and include the adsorption behaviour of gels, 
surface coatings, membranes in nano pores and biopolymers in restricted geometries. 

In particular, the formation of polymer networks at surfaces, is crucial for a number of 
applications, where surfaces have to be protected against forms of mechanical or chemical stress, 
such as abrasion and corrosion [^5]. Surface-attached networks also play an important role 
in several biomedical concepts, for example, to provide biocompatible, but stable coatings on 
implant surfaces Most theoretical treatments of these types of systems have been on the 
level of scaling theory, and analytical treatments have been lacking. Recently, Allegra and 
Raos investigated a confined polymer network, but modelled the effect of confining walls on 
a network by an harmonic potential and completely ignored the possibility of wall attachments. 

The aim of this thesis is to gain a better understanding, of a polymer network (or gel) that 
has formed at a confining surface. In particular, we investigate the simplest case where the 
confining geometry is two parallel planar surfaces. However, before we embark on the statistical 
physics theory treatment of this problem, some relevant polymer background should be given. 

1.1 What is a polymer? 

A polymer is a substance 2 composed of macromolecules, which have long sequences of one or 
more types of atoms or groups of atoms linked to each other by primary, usually covalent, bonds. 

1 Polymers in confined geometries is a multidisciplinary field, attracting research interest from scientists in 
chemistry, material science, biology and experimental, simulation, and theoretical physics. In 1998 the European 
Associated Laboratory (of the Institut Charles Sadron, Strasbourg and the Max Planck Institute for Polymer 
Research, Mainz) was created for the purpose of studying polymers in confined media. 

2 Although the terms polymer and macromolecule are used interchangeably, the former strictly refers to any type 
of polymeric material (rubbers, biopolymers, fibres, glassy and crystalline polymers) of which the macromolecule 
is the essential, common building block |47) . 
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Macromolecules are formed by the process of polymerization, that is, by linking together many, 
say N, monomer units through chemical reactions. The long chain nature of macromolecules 
is responsible for the characteristic properties of polymers and sets them apart from other 
materials. 

1.1.1 Flexibility 

Flexible polymers have a large number of internal degrees of freedom. The typical primary 
structure of macromolecules is a linear chain (the backbone) of atoms connected by chemical 
bonds and some pendant atoms or groups to satisfy the remaining valencies. By rotation about 
the single bonds in the backbone the molecule changes its shape or conformation. Since there 
are many of these bonds, a wide spectrum of conformations is available to a macromolecule 3 . 
The rotation of chemical bonds may be hindered by bulky pendant groups, so that some of the 
conformations become unfavourable. Sometimes the interaction between neighbouring groups 
leads to preferred sequences of bond orientations, which emerge as helical or folded sections in 
the molecules. Thus, a polymer is termed flexible if thermal motion is strong compared to the 
energy barriers associated with backbone rotation. 

1.1.2 Ideality 

The simplest measure of the length of a polymer chain is the contourlength L. This is the length 
of the stretched-out molecule, that is, for a chain of N bonds of length i the contourlength is 
N£. However, this length does not give a realistic measure of the size of the polymer chain, 
which in a molten state or in a dilute solution is coiled up. 

The conformational properties of long flexible chains can be described by the universal ran- 
dom walk model first introduced by Kuhn jSHJ- A chain is considered as a sequence of N 
randomly orientated bonds, each of length £. If the bonds are completely independent of each 
other, the conformation of the polymer chain resembles the trajectory of a diffusing particle 
under the action of a random force, for which the solution is well known |211 153] . If R is the 
end-to-end vector of the linear macromolecule, the mean square displacement is given by 



3 The degree of polymerization, N is usually more than a hundred. For example, polymerization of N ~ 10 4 
ethylene (CH2 = CH2) monomers will yield a giant polyethylene £j— CH2— \rf) macromolecule. Imagine now that 
this molecule only had three possible bond rotations; then the total number of shapes it may assume will be ~ 3 . 
Furthermore, the conformation is continuously changing due to thermal motion. A DNA molecule has N ~ 10 9 
links. Detailed analysis of these configurations is futile! In order to investigate the properties of a system of 
molecules, it is beneficial to consider a macroscopic system. In the case of macromolecules, even a single molecule 
is a macroscopic system with infinitely many possible conformations. Consequently, we can calculate relevant 
thermodynamic quantities of even single polymer chains by means of statistical mechanics. 

4 For non-linear, branched or star-polymers the radius of gyration R g , which is the root mean square distance 
of the segments from the centre of mass, is the approriate quantity, and is given by E? g = ^£ 2 N. 





Therefore, the characteristic size' 



i 4 R of the polymer is proportional to 2Va. The probability 
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distribution of an ideal chain endpoint to be at a distance R from the initial point is given by 
the Gaussian probability function P(R,N) oc exp (~ 2 3 jvi a ) [ a ^ so see Section l2~T] . 

For ideal chains, the finite volume of the segments and solvency effects are completely ignored. 
In reality, segments cannot overlap — called volume exclusion — which leads to chain expansion. 
The statistics of real, non-intersecting chains are described by self-avoiding walks instead of 
random walks 8 . 

However, polymers can adopt ideal dimensions in solutions in a so-called 0-solvent. In a 
good solvent, a chain expands from its unperturbed, ideal dimensions to maximize the number 
of segment-solvent contacts and the coil is said to be swollen. In a poor solvent, the chains will 
contract to minimize interactions with the solvent. However, competing with this effect is the 
tendency for chains to expand to reduce unfavourable segment- segment interactions (which is 
the excluded volume effect). If these two interactions are in balance, the polymer molecule will 
adopt unperturbed dimensions, and the solvent is said to be a G-solvent 8 . In short, if the 
concentration of molecules is high enough to be classified as a dense melt, the molecules will be 
forced to interpenetrate, so that (ideal) screened excluded-volume statistics may be assumed. 

1.1.3 The Entropic Spring 

The entropy of a macromolecule is described by the formula 

S = k B lnn, (1.2) 

where k B is the Boltzmann constant and SI is the number of possible conformations. When the 
chain is deformed, the entropy change AS is given by, AS = k B In fi^/fio, where the subscripts 
and A refer to the initial condition of no strain and the condition of the oriented state under 
stress, respectively. 

Near an impenetrable surface, the geometric restriction leads to a lower conformational 
entropy of the polymer. 

1.1.4 Elasticity 

Rubberlike elasticity is the consequence of molecular arrangements — in other words — chain 
flexibility. This is quite different from the elasticity of ordinary solids such as pure metals and 
crystals, where the elasticity or the resistance to deformation under external force, arises from 
the distortion of the intermolecular potential fields. For macromolecules under strain, the inter- 
molecular potential energy U remains nearly constant with or without strain [48) . The elastic 
driving force is therefore entirely from the tendency for macromolecules to randomize in order 
to attain the maximum entropy, and minimum free energy. For this reason, a macromolecule is 
often called the entropic spring. The entropy spring becomes stiffer at higher temperature T, 
since the tendency to randomize becomes stronger with more vigorous segmental Brownian 
motions like those of molecules in a liquid. This is in contrast with the behaviour of most crys- 
talline solids in which the potential energy is weakened and the stiffness is diminished at higher 
temperatures as a consequence of thermal expansion, which moves atoms further apart. 
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The above discussion on polymer elasticity can be summarized in terms of a thermodynamic 
equation of state for the stress /, given by 

^ \ de^j T \ cte ) T ( cfe ) T ' ^ ^ 

so 

where F is the Helmholtz free energy and e is the length of the strained sample. For an ordinary 
solid, like a diamond, the reverse of (Jd is true, that is, T (^) T is zero. 

1.1.5 Polymer Gels 

Gels are macroscopic network polymers that have three-dimensional structures in which each 
chain is connected to all the others by a sequence of junction points, called crosslinks. 

The crosslinks, together with the condition of chain flexibility, are responsible for rubber 
elasticity 5 . The classical theory of rubber elasticity j3Hj predicts that the free energy, T of 
deformation is proportional to the sum of squares of the principal extension ratios: 

i=x,y,z 

where N c /V is the crosslink density of the network sample. A major discrepancy of this model 
is that it assumes an affine deformation: if a macroscopic rubber sample is deformed by A, then 
the end-to-end vector R of any subchain between two junction points will be equal to A • R, 
after deformation. The affinity assumption implies that the crosslinks are spatially fixed, and do 
not fluctuate. In the James and Guth model the crosslinks are essentially unrestricted, and 
the resultant free energy in Q1.4JI is altered by a factor of ^. In 1975 a pioneering network model 
was introduced by Deam and Edwards [S], which models the effect of the network on a given 
chain by a harmonic localizing potential. In order to calculate the free energy of deformation, 
they resorted to the replica method from spin glass theory. The resultant Deam and Edwards 
free energy Q4.66JI . again shows the same strain-dependency than l|1.4|) . but with a different front 
factor and more terms depending on the crosslink density. In these phantom models, excluded 
volume interactions are ignored, and the polymer chains can pass through each other. 

In reality, polymers displaying rubber elastic behaviour, will deviate from phantom models 
like (|1.4j) mainly due to the presence of chain entanglements [3] . A detailed review of the role of 
entanglements and attempts to understand these topological constraints, like the slip-link and 
tube models, can be found in [T51 152] . 

Although network models are numerous 6 , the demand for theoretical treatments of surface- 
attached, confined networks, has remained unfulfilled. 

5 This is the simplest view of a material that is expected to contain network inhomogeneities |38| and other 
defects, like trapped entanglements, which will alter its elastic behaviour. 

6 A recent review, and experimental comparison between different models of unconfined networks is given by 

ESI- 



Chapter 1 — Introduction 



5 



1.2 Outline of chapters 

The thesis is structured as follows. In Chaper 1 we present a single, free macromolecule that 
is subsequently confined in box and plate systems. Within this simple single chain-system we 
introduce the concept of a stitch and deterministic bulk and wall crosslinkages. 

Chapter 2 deals with quenched disorder, due to random crosslinking, in the context of con- 
finement and a simple stitch-network model. 

In Chapter 3 the basic statistical crosslinking model for a confined, surface-attached model 
is presented. We extend the Deam and Edwards [0] phantom model, to include two parallel 
confining walls and the possibility of surface attachments or wall-links. This model treats the 
localizing effect of the crosslinks on the system, by a harmonic potential with a strength that is 
isotropic and strain- independent. In Chapter 4 we improve upon this assumption, by employing 
an inhomogeneous localization potential. For simplicity, this is done in the framework of a brush- 
network. In both models, we start by formulating the problem in the language of statistical 
mechanics and constructing the partition function. Thereafter, the free energy of deformation 
and the stress-strain equation are calculated. 

In this thesis we adopt the philosophy of first finding a tractable theory, under reasonable as- 
sumptions, before attempting more realistic situations. The first simplification, used throughout 
this work, is the assumption of ideal, Gaussian chains. The omission of excluded-volume effects 
(in contrast with the neglect of trapped entanglements) from the theory, is still a laboratory- 
attainable assumption to make. These flexible, intersecting chains are illustrated throughout 
the text by kinky, spaghetti-like lines, emphasizing the fact that we are dealing with theoretic, 
simplified entities. The second assumption, is that of a sufficiently crosslinked network. This 
underlines the fact that we do not work near the so-called sol-gel (second-order) phase transition 
region |24j . 

In the last chapter we briefly conclude the work and give an outline of possible improvements 
and future developments. 




Figure 1.1: (A) A schematic description of the preparation of a surface- attached network by 
simultaneous photochemical crosslinking and surface attachment. (B) The photochemical reaction 
used for surface attachment. [Picture taken from |4U| ,] 



Chapter 2 

The Single Confined Chain 



Before we can probe the physics of a polymer network formed at a surface, we have to formulate 
the problem in terms of a specific mathematical model. In this chapter we present the funda- 
mental concepts that we shall employ to construct such a model. A vital approach throughout 
is to use suitable Green's functions to investigate the statistical properties of the system. Using 
the Green's function approach, we firstly introduce a free Gaussian chain 1 , and then place it 
in a confining environment. The behaviour of a macromolecule is determined by the number 
of different conformations it can take. By placing the polymer in a restricted geometry, only 
certain specific conformations are selected. After confining the chain, we link its two endpoints 
onto the confining surface. These surface links define the simplest case of wall-links, and restrict 
the confined polymer to an even greater extent. Lastly, we introduce a single polymer-polymer 
link, a so-called bulk-link, under very specific conditions. 

The work in this chapter together with the prefatory treatment of disorder in Chapter |3] will 
serve as a platform for the full calculations. 

2.1 A free Gaussian chain 

One way of representing an ideal macromolecule is via the standard Gaussian or bead-spring 
model, Figure [23(b). For the standard Gaussian model, the chain conformation is specified by 
the set S = {Ro, Ri, • • • , R^v} of N beads, which can be thought of as N repeating monomer 
units of the chain. The conformations of an ideal macromolecule coincide with the random walk 
path of a Brownian particle. Since Brownian motion is a Markovian process, the ideal chain 
also belongs to the class of Markov chains. Let '(/'(Rn., Rn-i) represent the linear memory that 
describes the bonds between a pair of link neighbours. The memory of the chain direction is lost 
over a distance comparable to its persistence length 2 ~ £. The probability of a given polymer 



x We do not attempt to give a complete introduction to the mathematics of ideal chains, random flight polymers 
and the random walk analogy, as this will (and already does) expend many polymer theory textbooks |16M53l l£lll3 ) . 
A comprehensive study of the Brownian chain and associated probability law is given in |12| . 

2 The persistence length of a long, flexible, ideal chain is related to the effective Kuhn segment length £ of a 
freely jointed chain oi N = L/£ segments, such that (R 2 ) = Nl 2 . In the context of the Gaussian chain model, 
the Kuhn step length is defined by {(R„ — R n _i) 2 } = I 2 , and denoted by the term (RMS) link length. 
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(a) (b) 

Figure 2.1: (a) An ideal, long flexible phantom chain, and (b) its Gaussian model counter- 
part. The shape of the polymer can be represented by the set of position vectors of the beads 
S = {Ro,Ri, . . . ,Rjv}, or by the set of bond vectors {ri,r2, . . . ,rjv} ; where r n = R n — R n _i. 
The end-to-end vector R = Rat — Rq characterizes the size of the chain. 



shape S is thus given by VPfS 1 ] = n»j=i VKR-n, Rn-i), and represents the connectivity of the ideal 
macromolecule. For the bead-spring model, the bond lengths have the Gaussian distribution, 
such that ^(R n ,R n _i) is defined as follows: 

V(Rn,R„-l)= (^2) 2 exp j-J^Rn-Rn-l) 2 } • (2.1) 

The partition function is computed by integrating the distribution ^[S] over all possible con- 
formations [8]. 

When the chain is very long, the monomer index n may be regarded as a continuous variable. 
In this continuous representation, the discrete bond vector r n = R^ — R n _i is replaced by the 
functional derivative dH n /dn ^Bj. If the endpoints of the macromolecule are fixed, at say 
positions r and r' in space, the partition function is the Wiener path integral 

" R(7V)=r ™ n r 3 [ N (dR. n \ 2 

'R(0)=r' 

which is a Green's function and a solution of the following diffusion type equation: 



/•R(JV)=r , o fN / a-o \ 2 

G(ry,N) = Ml [m]exp{-^y o (-^) *,}, (2.2) 



\M-ji} Giiy:N)=S(r - T ' )SiN) - (2 - 3) 
The solution of (|2.2j) or (|2.3j) is the Gaussian distribution, given by 

G (^^» = (^) 3/2exp {-^^}- (2 ' 4) 

If all interactions are neglected, a free polymer chain is characterized by an average end-to-end 
length of R = (R 2 ) = £N? . The number of configurations of a free, ideal chain between points r' 
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and r and chain contourlength L between those points thus corresponds to the statistical weight 
of a random walk starting at position r' and ending at r, in N = L/£ steps, where I is the Kuhn 
step-length of a polymer chain. 



2.2 A Single Chain in a Box 

Consider a phantom molecule confined in a box of volume V = h 2 h z , with length h in the x 
and y directions, and h z in the z-direction. The statistical weight of the confined polymer chain 
which starts at r' and ends at r in N steps is given by the Green's function G(r ,r' ;N), written 
as a Wiener path integral JS] : 



f-H(N)=r o fN /qtd \2 

GMN) = Ml [P box R]exp{-^y o d n}, (2.5) 



l(0)= 

where the normalization N refers to the number of configurations, Q2.4JI of a completely free 
polymer chain, starting at r' and ending at r in N steps, with Kuhn step length £. The notation 
means evaluating the path integral only in the allowable box-region. In terms of an at- 
tractive potential A, the Green's function ()2.5|) can be written in terms of unconfined integration 

cR(JV)=r r 3 j-N pN 

/R(0)=r' 



G(r,r';N) = M [DR] exp i-— / B? n dn+ Ai^dn], (2.6) 

jR(0)=r' 1 2^ Jo J J 



where 

, , f 1 if R in V , 
A(R,n) = ln0(R n ) , with 6(R) = <^ (2.7) 

I otherwise 

represents an infinite potential wall. The solution of (|2.5[) is equivalent to solving the inhomo- 
geneous differential equation 

M ~ If S " A(r > N) ] G(r ' r ' ; iV) = " r ' )5{N) > (2 ' 8) 

which, with the substitution of the potential A in 1)2. 7J) . simplifies to a diffusion equation of the 
form given by equation ()2..3)) . together with the boundary condition that G(r, r'; N) = at the 
confining surface. The solution G is constructed from eigenfunctions that vanish on and outside 
the boundaries of the box 3 . The coordinates are separable 

G(r, r'; N) = G x (x, x'; N) G y (y, y'; N) G z (z, z'; N) , (2.9) 



3 The solution was obtained by using an eigenfunction expansion method |35|. but can also be obtained by 
using the method of images pj. 
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and since the polymer in a box is a symmetric problem, all the parts, for example G x , have the 
same form: 



G x (x,x';N) = -J2 



2 x— v . Px^X . PxKX 1 



sin • 



Px = l 



h 



sin ■ 



h 



02 ^2 



6XP ' ~ ~6~h? PxN 



(2.10) 



Equation (|2,9|) is the probability that one endpoint (n = N) of a box-confined polymer is 
r, provided that the other endpoint (n = 0) is at r'. The partition function of all possible 
conformations, for a polymer with endpoints free, is given by 



Z = J dr J dr' G(r,r';N) 



E 

Px,Py,Pz = l,3,. 



\h 2 h 7 



> ) ) / CX P 

Pip z y pi 



6 L h? hi 



(2.11) 



(2.12) 



The endpoints of the confined molecule in (|2.12j) are not fixed at positions somewhere in the box 
or at the boundaries, as will be the case in the next section. 



2.3 A polymer stitch 

Next we add two wall-links to the system by linking the ends of the polymer to the top and 
bottom planes of the box, creating a polymer stitch, as illustrated in Figure l2~2l We are interested 




f ^ ~ 6 ) C^j \ 



h, — e 



IH ft e ) 

V 2 ' 2' > 




Z = 



Figure 2.2: A single polymer chain "stitch" with end-points held at two fixed points 



in the partition function of a single polymer chain stitch in a box with volume V = h?h z . To 
this end, we fix the end points of the chain a small distance e from the floor and lid of the box 4 , 
permanently, at 

, ,h h . , ,h h , . , 

r' = (-,-,e) and r = (-, -,h z - e), (2.13) 



4 Instead of a polymer that is cross-linked exactly on the wall, we have localized the wall cross-link a distance 
e, of the order of the link length £, away from the walls, otherwise the partition function for the polymer vanishes 
as should be expected. 
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By substituting these coordinates in equation 1)2.9}) and 1)2.10)) , we obtain the partition function 
for a stitch confined in a box: 

Z(rJ;N) = Z x>y (N,h) Z z (N,h z ), (2.14) 

with its three Cartesian parts written as follows: 

/ 2 \ 2 ~ r £ 2 it 2 ^ 

Z x , y (N,h) = (-) Yl ex P{-65?^ + *4 (2 ' 15) 

Px,Py=1,3y 

Z z (N,h z ) = A £ sin(^) S1 n(^%^l)e--S P ^. (2.16) 

P*=l 

It is at this point necessary to look at the relationships that exist between the box dimensions, 
h z and h and the length scale of the polymer y/~N£. 

2.3.1 Three limiting cases 

The relationship between the confining box dimension, h and h z , and the size of the chain y/~N£ 
determines the type of physical situation we are facing. It will also dictate in which of these 
limits future calculations will be done. There exist three distinct limits: 



(i) VN£ < h g < h (2.17) 

(ii) h z < \^N£ < h (2.18) 

(iii) hg <C VN£ and h < VN£ , (2.19) 

depicted graphically in Figure l2~31 In the first case (i) ()2.17j) . we have a small polymer in a large 




Figure 2.3: Relationship between the size of the polymer and confinement dimensions: (a) A 
small polymer in a big box, 1)2.17)) : (b) Polymer big enough to be a stitch, 1)2.18)) : (c) A large 
polymer, densely occupying the box, ()2.19j) . 



confining environment, which intuitively means that it will be very difficult to fix the endpoints 
a distance h z from each other. This difficulty is portrayed by the sum argument in ()2.16)) . which 
in this case oscillates rapidly around zero. Physically, (i) is analogous to a short random walk 
in a large cavity |26j . Since the polymer is too small to fix at the top and bottom faces of the 
confining box, we investigate the original partition function 1)2.11)) . for a polymer with free ends. 
Since the characteristic chain size R is small relative to the box dimensions, the sums in the 
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partition function (|2.12j) will be dominated by ^ terms [Ej- I n this limit 5 , the free energy F is 
like that of a perfect gas, with corresponding pressure P, given by: 

dF 

F ^ -k B T In V, and P= — ^PV = k B T. (2.20) 

oV 

In the second relation (ii) 1)2.18)) the stretched out chain is small compared with h, but long 
enough to reach the endpoint cross-link locations in the z direction. The sum in Z(N,h z ) in 
(|2.16|) is then entirely dominated by the p z = 1 term, so that the /^-dependent part is given by: 



Z z (N,h z ,X) = l z Sin[^-} sin f^ e) ] e~^ N . (2.21) 

On the other hand, a continuous Gaussian approximation compares well with the sum over p x 
and p y in (|2.15f) . that is: 

/1\ 2 /2\ 2 7 7 I Nf2 ^iP 2 x+P 2 y) ^ 

Z x , y (N,h) * (-) (-) J J e { ^ } d Px dp y (2.22) 

o o 

= 2^N- (2 ' 23) 

which is the same as not restricting the x and y dimensions at all, and describing each of these 
components of the chain by a free chain Green's function (|2.4|) . In this limit of sufficiently large 
h, and infinitesimal e, the problem is equivalent to a chain restricted between two parallel plane 
surfaces with an h-independent partition function, 

z{hz) ~ siw^f exp V ~~6hf) ^ limit r z ^ ^ l ^ h \ ■ (2 ' 24) 

In the third relation (iii) (|2.19|) the polymer contourlength is large compared with both h and 
h z . Consequently, the dominating contribution to the sums in 1)2.1 5 J) and 1)2.16)) , is given by the 
Px = Py = Pz = 1 term, leading to a partition function (in the limit of small e) , 

Z{h,h z )~^jp ex p{-^V L (^ + ^ limit ( in ) [h s <h<VNi] , (2.25) 

which is dependent on all the parameters of the confinement. In this instance the single polymer 
chain will tend to fill up all the available space. 



2.4 Stitch Strains 

The mechanical properties of a polymer stitch may be investigated by looking at what happens 
during strain. It is possible to characterize the stress acting at a point on a substance in terms of 



5 If \/N£ < hi, the partition function Z ~ h' 2 h z = V, since exp{- \n 2 pj ( ) 2 } ~ 1, i = {x,y,z}. The 
remaining sum is then of the form ^Z^Li 3 l/p 2 > which converges to 



Chapter 2 — The Single Confined Chain 



12 



three principal stresses or extension ratios, acting along mutually perpendicular principal axes 
|48| I54j . The extension ratio Aj is defined as the deformed length in direction i divided by the 
original length of the sample |47| . 

In principle there are many possible deformations that can be employed to verify the pre- 
dictive ability of theoretical elastomer models describing stress-strain behaviour. However, in 
practice it is uniaxial deformation that is most often studied, due to its experimental simplicity 6 . 

An important characteristic of macromolecules in the gel state is that any type of applied 
stress will most readily influence the gel's shape, without appreciably changing its volume 
|4fi| I48j . This observation has led to the terminology incompressible or indilatable to describe the 
mechanical behaviour of gels. In most cases we shall thus be concerned with an isovolumetric, 
uniaxial deformation (|2.2fij) . where an elongation (compression) by a factor A along the z axis 
results in compression (dilation) by a factor l/\/A along the other axes. The tensor A expresses 
the isovolumetric, uniaxial (macroscopic) deformation of a point r — > r' = A • r, 

/ A-2 \ 

A = A-5 , A > 0. (2.26) 
V OA/ 

The fundamental property of polymer gels, namely the constancy of volume during deformation, 
makes it possible to define the complete state of strain in terms of a single parameter A. 

2.4.1 Case 1: Deforming the plate system 

Firstly, we investigate the effect of strain, defined by A, on a macromolecule linked at fixed 
positions to two parallel plates fFigure l2~4"|) . This situation corresponds to limit (ii) with h — > oo. 
In this case it is appropriate to use the partition function Z, already derived in (|2.24l) . to calculate 




A 



Aft. 



i _ ih h \ 

V2' 2' > ' / — ( h h A / 

r - U.a> 9,/i> ( ! 



(a) (b) 

Figure 2.4: Polymer stitch in the (a) the unstrained state; (b) strained (unidirectional com- 
pression) state, confined between two long parallel plates, corresponding to limit (ii). 



6 However, only general 6i-axial strains cover all accessible pure homogeneous deformations for an incompress- 
ible material, and is the preferred method of testing real network theories that account for topological features 
like trapped entanglements |49j . 
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Figure 2.5: PZoi of the stress f (|2.31j) versus the deformation ratio A /or a p/aie system, in 
limit (ii). The solid line is for a stitch confined between the plates; the dashed line corresponds 
to a free chain (with fixed endpoints). Since the free energy was evaluated in the limit (ii), only 
small elongations, Xh z <C y/~Nl are allowed. 



the free energy of the plate-system as follows: 

F = —k B T In Z(hz) (2.27) 
= —k B T In 



Tre 2 1 k B TNit 2 l 2 



3N£ 2 B 



+ ^^. (2.28) 



If we perform a simple elongation (A > 1) or compression (A < 1), the change in free energy is 
given by: 

AF _ wlnA+ *™p (2 . 29) 

which has similar characteristics to that of De Gennes' scaling theory calculation |Sj for the 
energy required to squeeze an ideal chain (of unperturbed size Rq) trapped in a tube or cavity 
of diameter h z : 

R 2 

A F deGennes = k B T-^ ,Rq = Ni i. (2.30) 

The free energy of the form in (|2.30f) holds generally for any type of ideal chain confinement |26j . 
It originates from the decrease in the number of available conformations — or lowering of the 
entropy — when the polymer is restricted. The only force acting on the system during strain is 
the tensile force in the direction of extension (or compression) A. The magnitude of the force 
per unit cross-sectional area measured in the unstrained state, depicted in Figure 12*31 is given 
by 



V dX V \X 3h 2 A 3 



An alternative way to model a stitch is to simulate the effect of the wall confinement by a simple 

harmonic potential [ITj. If the potential is U = H 2 dn, the wall-confinement is enforced 

_ i 

by choosing the localization parameter to be of the form q ~ h z 2 . The partition function in 
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this case is analogous to the propagator 7 of a particle in an harmonic potential, with initial 
and terminal points of its trajectory fixed on one of the plates |19| . The elastic free energy 
calculated from the harmonic approximation is consistent with (|2.29|) . and confirms the form of 
the stress-strain graph, Figure 12,51 

When an unconfined ideal chain (|2.4j) is deformed from an end-to-end length h z to Xh z , 
its end links will undergo an external stretching force of mag nitude / = k B T ^ JiUHl- The 
important difference between a confined and unconfined polymer stitch, is illustrated by the 
stress-strain relationship in Figure 12*31 

2.4.2 Case 2: Deforming the box system 

Secondly, we investigate the effect of strain, defined by A, on a macromolecule in the limit (iii), 
y/~Nt 3> h, as shown in Figure |2~E1 In this limit it is appropriate to perform calculations using 
the partition function already derived in (|2.25|) . 




It is found that the work done to strain the system, 

AF = 5 l '™' V (¥ + (^4)' (2 ' 32) 

is both h and h z dependent. The only force acting on the system during strain is the tensile 
force in the direction of extension (or compression) A. The stress-strain relationship is given by 

_ 1 d(AF) _ k B TN£ 2 Tr 2 ( 1 1 \ 



and is only valid for small elongations Xhi -C yNl. Both elongation (A > 1) and compression 
(A < 1) is represented by a single curve on the next page (Figure l2~7]) . illustrating the 
theoretical relationship between force (stress) and A when relation (iii) holds. 



7 Similar calculations that employ a Green's function solution for a harmonic localization of the crosslinks, is 
given in more detail in Chapter |I] 
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/ 

1 Elongation 




Figure 2.7: Plot of the stress (f) (I2.33|) versus the deformation ratio (X), in limit (Hi) when 
h <C y/~N£, for a long polymer confined in a box. 



2.5 One bulk-link 

Consider two, identical macromolecule stitches, each of length L = N£, confined between two 
parallel plates, in Figure l2~HT a). Next, the chains are joined together permanently at position R, 




(a) (b) 

Figure 2.8: Two polymer stitches, (a) without a junction point; (b) with a junction point 

at R = (X, Y, Z). The wall links are situated at permanent points, v\ = [x\,y\,h z — e), 
r2 = (X2,y2,e), r 3 = (x 3 ,y 3 ,h z - e) and r 4 = (x 4 ,y 4 ,e). 



such that both chains meet at half of their respective contour lengths, |iV. We obtain a system 
consisting of four wall-links and one bulk-link, as shown in FigureEBl(b). The partition function 
for the system is thus the statistical weight of one large star polymer 8 of total length 2L. The 
statistical weight is determined by averaging over all possible positions of the junction point R 
for the four half-chains JS], that is: 

Z = J dRG(r 1 ,R^)G(R,r 2 ;y)G(r3,R;y)G(R,r 4 ;y), (2.34) 

8 A regular n-star polymer is a macromolecule containing a single branch point from which n linear chains (iden- 
tical with respect to constitution and degree of polymerization) emanate |26) . The size of such a macromolecule 
is described by the mean square of the radius of gyration: -R 9 4 arms = § N£ 2 /6. 
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with the points ri, r-2, r3 and r-4 at fixed positions on the plates, shown in Figure l2~%l (b). 
The Green's function G(ri,R; ^) represents the statistical weight of a chain portion which 
starts at ri and ends at R in y steps. Since the chains are only confined in the z dimension, 
each Green's function G is the product of two free chain contributions, (|2.4|) and one confined 
statistical weight contribution (|2.1U|) . for example: 



N 

G(ri,R; y ) 



N N N 

G x (x u X; —) G y (y u Y; -) G z (z u Z; -) 



ir£ 2 N 
2 



cxp 



Nl 2 



((X- Xl f + (Y- yi ) 2 ) 



EirpZ ftp(h z 
sin — — sin — 
,, z h z h z 
p=i 



N 



(2.35) 



(2.36) 



If this one-link system is strained by simple elongation A, as defined in (|2.26j) . the change in 
free energy is given by 



AF 



-k K T In ■ 



Zi 



(2.37) 



where Z\ denotes the strained and Z\ the unstrained states of the system, defined in (|2.34l) . In 
the limit (ii), y/~N£ > h z S> e, equation (|2.37f) reduces to: 



AF ~ k B T < 71n A + 



N£ 2 



£tt 2 N / 1 



3h 2 



x 2 ~ l 



+ 



i=l 



- - 1 ) ) (2.38) 



with d 2 depending on the chosen, albeit fixed, relationship of the wall linkages. The constant 
d 2 is translation invariant and ought not play a dominant role, since the problem is symmetric 
in the x and y coordinates 9 . However, for the above expression to be valid in the limit (|2.18|) . 
one must limit the distance between wall links for this two-stitch system. This can most easily 
be seen, when the number of coordinate constants xi and yi is reduced, such that x\ = X2 = x a 
and X3 = X4 = Xb- In this case d 2 = (x a — Xb) 2 + (y a — Ub) 2 , and gives a measure of the distance 
between two wall-linkages on each plate. The entropy of the system will thus increase as the 
distance between points r a and r& decreases (Figure l2.9|) . Conversely, when d, or the distance 
between wall-links in the x-y plane, increases, the chain portions are more stretched out and 
fewer configurations are available to the polymer. 

When the stress / is calculated in the manner of (|2.31j) , we obtain a stress-strain relationship 



9 The deformation factor (1/A — 1) after the constant d 2 is zero if the plates are simply stretched in the z 
coordinate by a factor A. In Equation 12.3811 d 2 is present since an (optional) isovolumetric deformation was 
enforced. 
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Figure 2.9: Two polymer stitches with a junction point atTl= (X, Y, Z). The top and bottom 
wall links are situated at permanent positions with equal x and y coordinates, such that the left 
links have coordinates (x a ,y a ,e) and (x a ,y a ,h z — e), and the right has coordinates (x&,yh,e) o,nd 
(xt,,yb,h z — e), for top and bottom plates, respectively. The arrows show movement of the fixed 
links to a common position that correspond to maximum entropy. 



similar to (12.311): 



/link 



k B T 
V 



2Ntt 



2e2 



i 



3d 1 



Nt? A 2 



(2.39) 



The coefficient of the confinement term in (|2.39j) , 2N ^ l -/ has increased by a factor two compared 
to (|2.3H) . which portrays the relative difficulty in deforming a system consisting of two chains 
with an added bulk linkage. Also, if we let d 2 = 0, we recover the case of only two wall-links. 
Since d 2 -C v 7 ~Nl, the second term in (|2.39j) will dominate, such that the stress /ii n k will not differ 
much from / for the simple stitch system (Figure I2.5j) . This is illustrated by a similar stress- 
strain plot of (|2.39j) . shown in comparison with the classic stress-strain curve for an unconfined 
"network" of only one bulk-link (and no wall-links) |48| . 



/link 

2r 



-2 



Classic rubber elasticity 



2.5 



f— Confined single bulk-link 



Figure 2.10: Plot of the stress (TiinkJ (|2.39|) versus the deformation ratio (X) for the confined, 
single bulk-link plate system (solid line), in limit (ii). The dashed line is the stress-strain curve 
corresponding to the classic rubber elasticity model, with f oc A — 1/A 2 , for only one bulk link. 
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2.6 Contemplating stitches 

In this chapter we briefly introduced the basic concepts which play a part in the formation of 
a network at a surface: the Gaussian chain, confinement, homogeneous strain, wall-links and 
lastly a non-random, single bulk link. 

Confinement reduces the number of possible chain shapes, and decreases the entropy. In a 
confined chain, the role of the characteristic chain dimension Ro is reversed, relative to its free 
chain counterpart. This leads to a chain which does not obey the Hooke "law" for stretching, 
so typical of an ideal chain. 

Since the Green's function of confinement included intractable sums, we were forced to work 
in certain limits, most notably, h z -C yN£. Consequently the stitches were not allowed to 
stretch too much. This restraining condition also applies to an unconfined phantom chain. Its 
force-extension relation is subject to the same limitations as the Gaussian distribution function 
from which it is derived. For end-to-end extensions approaching the total contourlength of the 
chain, the Gaussian approximation b ecomes progressively inaccurate, and. will fail if A max ~ 

It might be a reasonable, first approximation to simulate the real network (discussed in 
the Introduction) by a network of stitches: two parallel planar surfaces stitched together, with 
the possibility of link-formation in the bulk region between the plates. If we envisage such a 
network, restricted between two walls, we should work in the limit (|2.18j) since the stitch-system 
has only confinement in one coordinate (as opposed to a box-confinement). Each chain must 
also be large enough so that it can be fixed to opposite walls with ease. These links are formed 
an infinitesimal distance e from the walls. It was shown, for example in Section f2.41 that e does 
not influence the macroscopic properties of the system. 

During network formation, crosslinkages are formed at random. However, in this chapter we 
placed them at completely determistic positions. This descrepancy will be attended to in the 
next chapter. 



Chapter 3 

The Stitch Network 



How does one model a system with permanent, but random, constraints? 

In this thesis we are investigating the properties of a polymer gel, formed at a surface. 
The gel consists of long macromolecules that form cross links with the surface, in addition to 
the polymer-polymer links. We specifically consider a gel confined between two parallel plane 
surfaces. This symmetric problem is related to two plates, which are stitched together by very 
long chain threads as shown in Figure mi In the previous chapter we isolated one long thread, 
and fixed its terminal ends to the top and bottom confining plates. This was called a stitch. 
We included a mock bulk link that had its exact location on the chains, chosen beforehand. In 
a real network, the precise location of the linkages is not known beforehand, and will tend to 
vary from one sample to the next. Any realistic network model should thus take account of this 
randomness. 

In the present chapter we shall apply the replica method to the stitch network model to show 
how one can describe and handle random cross-linking mathematically. 

3.1 A random bulk link 

The first crucial step is to rectify the assumption of a given bulk link (Section 12. 5j) by making 
it random. Let s measure the arc length, s £ [0, L], along a Gaussian chain of contour length 1 
L. We simplify the toy model of stitches by having specific, but arbitrarily placed wall links. 




Figure 3.1: A polymer network of stitches, confined between two plates of spacing h z . 

1 In terms of the notation previously used, s = n£, such that L = N£, where £ is the average link length 
introduced in Section f2. 21 
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They are always chosen to be at the terminal points, corresponding to s = and s = L, 
of the polymer chain. We choose to divide the stitch network into a macroscopic number of, 
say M, subsystems, each having one random bulk link. The total free energy of the stitch 
network system would then consist of the sum of the free energies of the M subsystems, plus 
a contribution that comes from the interactions of the subsystems. For a system of phantom 
chains, we shall ignore any interactions. 

We continue by isolating a subsystem, a portion of the stitch network, as follows. Consider 
two stitches, each of length L = N£, confined between two walls of spacing h z , illustrated in 
Figure E21 



! 




Figure 3.2: A portion of the stitch network, as seen in isolation, consisting of two stitches 
linked at spatial position R. 

The wall links are formed an infinitesimal distance e from the walls 2 , at spatial positions ri, 
T2, r3 and r4, which do not vary from sample to sample. The bulk link is formed at a random 
arc length location s a and Sb on the two chains. This random linking is the source of the disorder 
in the system. Intuitively, any relevant observable of the stitch-system should depend on some 
general averaged characteristics of the random crosslinking. 

The main goal of this chapter is to introduce general methods for dealing systematically with 
the statistical mechanics of a random system. For the simple case of one bulk link, it might not 
seem worthwhile to employ strategies like the replica method (Section I3.2.3JI . This is especially 
the case when the crosslinkages are distributed uniformly. 

In the sections that follow we adapt the basic Green's function approach, developed in 
Chapter 12 to a system with disorder. 



2 The parameter e does not influence the physical properties of the system. For example, the stress-strain 
relationship is entirely independent of e, as investigated in section l^7H 
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3.2 Quenched disorder 

In Chapter [2] we used standard statistical averaging to compute the free energy of the system at 
a given cross linking position. For a real network, the linkage positions {S} on the chains are 
unknown, but we will assume that they are variables which vary randomly from sample to sample 
- called quenched variables - distributed according to a probability distribution P(S). Statistical 
physics averages that are generally employed, should be altered for quenched variables. 

3.2.1 Green's functions 

The partition function, or statistical weight, of the system with its degrees of freedom frozen 
at a particular set S = {s a ,s b }, and with its junction point at a specific spatial position R, is 
given by 

Z (R, S) = G (n, R; s a ) G (R, r 2 ; L - s a ) G (r 3 , R; s b ) G (R, r 4 ; L - s b ) . (3.1) 

In Equation ()3.1|) we employed the composition property of the Green's function to divide 
the total contour length (= 2L) into four smaller chain intervals. This is done analogous to 
Section [2.51 but with one crucial difference, namely that the interval lengths are not pinned 
down at L/2. The Green's function (|2.36f) of one chain section, with unknown radius length s a , 
is thus given by 

G(ri,R;a a ) = G x (x u X; s a ) G y (y u Y; s a ) G z (z u Z; s a ) (3.2) 



'v 

3 f 3 

2ttIs^ 6XP 1 27 

oo 



£t[(x-«r+<r-*r]} 



2 ^ irpZ np(h z -e) -^s a . . 

x — > sin sin '- e . (3.3) 

h z ^ h z h z 
P =i 

In order to obtain the partition function of all possible conformations of the system with a 
particular set of quenched variables, we have to average over all possible spatial positions of the 
junction point R: 

Z(S) = J dRG(r 1 ,R;s a )G(R,v 2 ;L- Sa )G(v 3 ,K;s b )G(K,v i ;L-s b ), (3.4) 

It might be possible to compute the chain configuration, for a specific choice of the quenched 
variables s a and s b , that minimizes the free energy of a specific sample. However, in a general 
disordered system such a calculation would be deemed intractable, since the number of quenched 
variables per sample is usually much larger than two 3 . Luckily, methods of statistical mechanics 
often prove fruitful in the limit of large systems. 



3 For example, in the Edwards- Anderson model |15| for a spin glass the interaction Jij between every spin 
pair <Ti and Oj is random. For N spins there will be ^N(N — 1) quenched variables per sample. Even numerical 
averaging over different realizations is too time-consuming |34j . 
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3.2.2 Averaging 

In the statistical mechanics of a random system, there is a need to perform two distinct averages 
|2L)| . Firstly, there is the thermal average that is carried out for every possible conformation, in 
a specific realization of {S}. This kind of average was performed in l[3.4|) . by means of a Green's 
function approach. Secondly, there is the disorder or configuration average, which is an average 
over all possible distributions of the quenched variables. 

In order to make the scenario of the random linking more general, the thermal average may 
be rewritten as follows. Let Tis denote the generalized Edwards Hamiltonian of a network 
subsystem composed of two randomly linked stitches, 



n S {{IU}}/k s T = V4 [ L ~R-(si)d Si - f L A [Rj(sj)] ds , 
i=l 1 ■'O J 



(3.5) 



where Rj(s) is the position vector of the ith. chain segment at the arc length s. The Hamiltonian 
describes the connectivity and ^-confinement of the chain system |22j . Although the four chain 
pieces seem to emanate from a common junction point at R (Figure l3.2|) . the respective arc 
lengths are not independent of one another, since they were formed from just two stitches. This 
adds a constraint to (|3.5j) : si = s a <^=^ S2 = L — s a and S3 = Sb <^=^ S4 = L — Sb- For a 
specific sample, the partition function is expressed as the product of separable path integrals, 



r /■R-l(0)=r 1 /-R 2 (i-Sa)=R />R 3 (0)=r 2 /•R 4 (L-s0=R jL n s [Rj] 

Z s = dR / / / [JjDRifofle 

J Jb 1 ( i „)=R JR 2 (0)=r 2 JR 3 ( Si) )=R JR 4 (0)=r 4 i=1 



fc R T 



(3.6) 



which is equivalent to the previous Green's function expression (|3.4[) . The free energy of a system 
in some particular crosslinkage state S = {s a , Sb} is the logarithm of ()3.6|) : 

F (5) = —k B T In Z (5) . (3.7) 

Let the probability of finding this particular state be P5. Then for a canonical ensemble the 
effective free energy is given by 

^ = ^FsF(5), where ^P 5 = l. (3.8) 
{<S} {S} 

The effective free energy is the observable or experimental free energy. A correct theoretical 
calculation of the average value of T over the entire ensemble is expected to correspond to its 
experimental value. The sample free energy F(5) is an extensive quantity and known to be 
self-averaging [201133! • We can thus expect that variations in F(<S) from one sample to the next 
and deviations from the average value J-, will go to zero in the thermodynamic limit V — > 00. 
In order to obtain an experimentally relevant quantity, only variables with the self-averaging 
property should be disorder-averaged. 

Taking a sample average of the partition function (|3.6|) instead of the free energy, is known 
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as the so-called annealed case: 

•^annealed = ~k B T In [E{s } P S 2s] , (3.9) 

and would give the wrong physics, since the crosslinkage in ()3.9|) is allowed to change in response 
to the chain conformations. In reality, {5} is fixed for each sample. Even if the substance 
between the planar surface is deformed, the crosslinkage for a particular sample will stay fixed 
in response to the deformation. Beyond what has been said, the partition function Z$ is also 
not suitable for averaging since it is not self-averaging and not physically observable. 

The correct strategy is thus to average the self-averaging free energy F (S) over the distri- 
bution of the permanent, but randomly chosen arc lengths: 

F=[F(S)] S = E {s} PsF(S) (3.10) 

= -k B T f ds a [ ds b P s In Z(S). (3.11) 
Jo Jo 

The configuration average is denoted by [. . .} s . Since the free energy for a specific S agrees with 
its configuration average [-F(5)] s in the thermodynamic limit, it seems that either quantities 
may be used in further calculations. However, the average [-FX5)]g is more suitable to work 
with, since it is a physically measurable quantity 4 . 

3.2.3 Replicas to the rescue 

Performing the average (|IS.11(I is usually 5 not feasible, because the dependence of In Z on S is 
often very complicated. In order to address the averaging the following identity, 

lim 2 - 1 = InZ or ^-Z n \ n=0 = In Z (3.12) 
n->o n on 

was first employed by |15M44| in theoretical models to predict the spin glass phase in disordered 
magnetic systems. This technical trick (|3.12jl . called the replica method, is useful because it 
is often easier to evaluate [Z n ] than [In Z\. Consider the partition function (j3.fi j) taken to the 
integer power n: 



Z n (S) 



n / dn(a) n / vn 



6 a) 



ex P \-^2n s [R^]/k B T\ . (3.13) 

a=l 



The above quantity is the partition function of n independent identical copies or replicas of the 
original system. In ()3.13j) the subscript a labels the replicas, and not the quenched variables S. 



4 The free energy is not the only self- averaging quantity. For example, the neutron scattering form factor is 
also a useful measurable quantity to average |52| . 

5 An example of a model where the disorder average has been solved analytically, without implementing replicas, 
is the random-energy model |10llll| . In this model the independent energy levels {Ei} are the quenched variables, 
and the sample average of the free energy is calculated by a microcanonical argument. The infinite range Ising 
model with p-spin quenched random interactions is also exactly solvable, and the results confirm the replica 
method 1271. 
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The chain conformation may fluctuate from replica to replica, but the cross link constraints are 
fixed for each sample replica. The scheme of the replica method can be described as follows. 
Firstly, the quenched average of the replicated partition function (j3,13|) for integer n must be 
calculated. Next, the analytic continuation of the the resulting function, [Z n (<S)] (S , should be 
made for an arbitrary non-integer n. Lastly, the limit n — > should be taken. In the end, the 
effective free energy is uncovered: 



The replica method might initially seem mathematically dubious since one evaluates [Z n (<S)] 5 
for integer n and then extrapolates the result to n — * 0. In the pioneering decade of spin glass 
theory, the replica method was thought to be the cause of some controversial and inconsistent 
results. In particular, the free energy of the Sherrington and Kirkpatrick (SK) model |H] was 
shown to be unstable and exhibit a negative entropy in the low temperature (< T c ) region. 
However, these unphysical results were a consequence of the replica symmetry ansatz and possi- 
bly the incorrect reversal of limits [HI, and not due to the replica procedure itself. Alternative 
mathematical efforts, including replica symmetry breaking schemes and the mean-field TAP 
(Thouless, Anderson and Palmer) approach 6 , agree with the SK solution at and above the criti- 
cal temperature T c . In all cases where the calculations can be performed by a different method, 
the replica approach is confirmed and gives sensible results ^01 I2Z| • 

In polymer network theories the use of the replica method is simpler and undisputed [21 
I37| I36j . because there are no competing interactions ("frustration") or need for replica symmetry 
breaking as in the case of spin glasses 7 . 

6 Comprehensive literature and reviews exist on the subject of the theory of spin glasses 20 34 [3]. Recently 
it has been shown that the TAP and replica methods are equivalent in the SK model 0- 

7 An example of where various methods are used to probe the disorder in a polymer system, is the problem of a 
Gaussian chain trapped in a medium with randomly frozen obstacles |17|. Variational calculations are performed 
with and without the breaking of replica symmetry, and shown to be consistent. It is also shown that the annealed 
and quenched cases are two distinctly different situations. 

Furthermore, many network theory problems have been approached without the use of replicas |45| . and corrob- 
orate the replica results. An example is that of determining the neutron scattering function of polymer networks 




(3.14) 
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3.3 The actual calculation 



In this section we apply all the available tools — the Green's functions of Section f3 . 2 . 1 1 and the 
replica method — to a subsystem of the stitch network (Figure l3.2|) . The replicated partition 
function of the single bulk-link system can be abstractly expressed as: 



Z n 



da.® gv> 



(3.15) 



where the Green's functions in (|3,2j) are also replicated 
g (a) _ G (a) ( ri)R W ;Sa ) ( R W, r2 ;L- So ) G^ (r 3 ,K^;s b ) G^ (R^,r 4 ;L- Sfe ) . (3.16) 



The first step of the replica procedure involves the configurational average of Z§ (|3.13j) . which 
in terms of Green's functions, is given by 



ds n 



ds b Yl 



a=l 



P(s a ,S b ) . 



(3.17) 



For the sake of illustrating the replica method in the simplest manner, we choose P(S) to be the 
uniform distribution, that is, P(S) = 1/L 2 . Although the bulk links are uniformly distributed, 
we do not allow them to form exactly at the ends, s = and s = L, where the wall links are 
situated. Since the location of the wall links of a stitch network is by definition chosen to be 
non-random, the x and y coordinates of the stitch ends may be taken to coincide at the top 
and bottom plates respectively 8 . In the limit y/~N£ S> h z , the first term in the eigenfunction 
expansion ()3.2|) , corresponding to p^ = p% = p^ = p^ = 1 in (|3.16f) , dominates the partition 
function. In the case of the uniform distribution of the quenched variables, the configurational 
average of the nth power of the partition function is 



L 2 



(const) 11 



L-e 



cxp 



ds a 
3n L 



L-e 



ds b [s a (L - s a ) + s b (L - s b )] 

jxi-x 2 ) 2 + (yi-y 2 ) 2 ' 
s a (L - s a ) + s b (L - s b ). 



(3.18) 



where "const" refers to a constant factor, dependent on the wall spacing h z , but not on the 
disorder. The e-factor in the integration, mathematically prohibits the bulk-links from being 
formed at the walls, which would result in double crosslinking. 

3.3.1 The method of steepest descents 



It is convenient to transform to a new set of integration variables, a x = \ s a — ^ \ and a 2 = \s 
where Oi € ■§)• Let the measure of the distance between points on the top and bottom 



L 

2 I' 



plates be denoted by d 2 , analogous to Section l2~5l d 2 = (x± — x 2 ) 2 + (yi 



2/2 J 



After the 



This was also discussed in Section 12.51 and illustrated in Figure |2 
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transformation, the exponential in (|3,18|) only has terms quadratic in 

1 



L 2 



(const) 11 



da, / da 2 e s(^) , 



(3.19) 



with the exponent g defined as follows 



9 {<Ti,<t 2 ) = nln 



T 2 

^ 2 2 

a, — cr„ 

2 1 2 



+ 



3nLd 2 ( L 



2f 



2 2 

2 1 2 



(3.20) 



If the dominating factor 3n ^ d in g is large, it is possible to evaluate the integral by the method 
of steepest descents or saddle point approximation |28| . The point (a*, crj ) where g is minimized 
is found by making g stationary with respect to a 1 and a 2 ; 93 %l'° 2 ' > = and 9b 1 = 0. 
A global minimum is found at (cr*,a 2 ) = (0)0)) only if L < 3d 2 /£. In terms of the original 
arc length coordinates, the system will thus favour a bulk link to be formed at s a = L/2 and 
s b = L/2. 

Next, the function g (a 1 ,a 2 ) may be expanded in a Taylor series around (cr*, a 2 ): 



z)rr 2 



g(a l7 a 2 ) = g(a* 17 a*) + ±{(a a - a?) 2 + 2 (a, - a*) (a 2 - a*) 

] -} + ... (3.21) 



+(^-^) 2 8g tr 2) i' + 

3rad 2 , 1 o 2n /3<i 2 



+ nln -L + -p- ( — L ) (of + o" ; 



(3.22) 



We choose to work with very long macromolecules, such that the contour length L is very large 
(but finite) throughout the above Taylor expansion. In this limit, the higher order terms in 
(|3.21|) are insignificant, so that (|3.19|) becomes 



I 2 3nd 2 f f 2n I 3d 2 T \ / „2 i„2\ 

[Z n \ s = j ? (constp) n e-^ / d^ / da 2 e~^ ( — " L) 



1 / 2 3n d- 

£2 (const -g) e a 



do-! / do- 2 ^ ■ 
k=o 



k! 



2n /3d 



L (of + o- ; 



(3.23) 



(3.24) 



However, since the integration limits are finite, we are forced to revert to a series representation 
of the exponential in the integrand (|3.23j) . The average of the replicated partition function 
■24|) . is then given by 



1 — n 



'vr 2 L Ad 2 
+ 



3k 2 , 



L 



In 



6e^ 

W \2l) L 2 



+ 0(> n 2 ) terms (3.25) 



Lastly, by taking the replica limit n — > as in (|3.14|) . the effective free energy T is obtained 

(3.26) 



m £vr 2 L 4d 2 , 
f/fc B T ~ — — — + — - In 



3h 2 



L 



r6e 4 




) 4 - 


1 




\2£s 


' L2 - 


~3' 



which is only valid for large, but finite chain length L, such that \J £L 3> h z , and when the 
relation L <3d 2 /£ holds. 
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3.4 Deforming the stitch network 

In Section 13.11 we defined the free energy of the phantom stitch network as a whole, to be the 
sum of M subsystems. Each single-bulk link system has an effective free energy T 1)3. 26 j) . If a 
stitch network is fabricated from long, but finite length chains, and the relation h z <C L < 3d 2 / £ 
holds, the replica procedure results in the following elastic free energy: 



, _ U^ 2 L Ad 2 , 

•^network ^ k B T M < ^ + -jj- - In 



6e 4 /3\ 4 2 



hi \2£ L 2 



(3.27) 



The first contribution to the free energy, ^ is ascribed to the i-confinement, and is typical 
for an ideal chain trapped in a one dimensional cavity of diameter h z [22j. The second term, 



4*j-, corresponds to the free energy of four unconfined, ideal chains (in two dimensions) with an 
end-to-end distance d. The last term in 1)3.27)1 comes from the normalization factor. In the case 
of a uniform distribution, P (S) = 1/L 2 , it is possible to find a similar free energy, without the 
replica method, by minimizing F(S) (|3.7|) with respect to s a and Sb- 

Using the above free energy (|3.27|) . it is possible to investigate what happens during strain. 
Let the stitch network between plates, be strained by an isovolumetric, uniaxial deformation A 
(|2.26|) . Rj — ► A • Rj. The stress as a function of the deformation ratio A, can be calculated as 
in Section E" 



k B TM (7 2tt 2 £L 1 id 2 1_\ 

/network- y ^ 3h 2 X 3 £ L \ 2 J ' [ ' 

and is found to resemble the stress-strain relationship for a network of stitches that has an 
average number of M non-random bulk-links, chosen beforehand to be at (s a , s&) = (-§,-§)• 

The stitch network model served as an exercise in the application of the replica method, and 
as motivation for the search of a better network model. 



Chapter 4 

The Real Confined Network 



In this chapter we shall extend Deam and Edwards' pioneering model of network formation j^j, 
to include a confining surface and random polymer-polymer and polymer- wall links. 

A polymer network formed at a surface is a system governed by a variety of constraints. 
The surface at which the gel forms, plays the role of a confining geometry that restricts the 
network chain conformations and positions. During network formation, the macromolecules 
form permanent links with one another as well as with the confining surface. Since the links are 
unbreakable, the junction points will be localized to a certain extent dependent on the crosslink- 
density. In a real network the exact location of the links is unknown prior to fabrication and 
varies from one specimen to the next. However, any two samples having different cross-link 
realizations are expected to exhibit similar macroscopic properties. 

In previous chapters we investigated single-chain systems subject to similar conditions of 
confinement and random linking. This was done to provide the basic groundwork for under- 
standing and approaching the actual network (many-chain) problem. In the present chapter we 
traverse through the "making" of an ideal polymer network formed at a surface, starting from 
a group of chains to the finished gel-product 1 . 

4.1 The Collection of Phantom chains 

Prior to linking, we have a system of M independent chains. They are assumed to be long, 
linear and have no charge. The calculations in this chapter are performed on phantom (pure 
Gaussian) chains, thus neglecting any effects due to inter- and intra-molecular forces. The term 
"phantom" emphasizes the fact that the assumption of no steric effect is clearly unphysical. In 
reality each monomer segment of a chain will possess some volume which is excluded from other 
segments (HEH; anc ^ conse( l uen tly the polymer will swell. However, in a dense polymer solution 
the excluded volume interactions are screened, and each chain is essentially ideal and Gaussian 
[S]. The collection of phantom chains may thus be visualized as a rather dense, overlapping 

Note: this is done through the eyes of a theoretical physicist, not a chemist. In referring to an "ideal 
network" , we mean a 3-dimensional network of phantom chains, joined together and fixed to the walls by random, 
unbreakable chemical bonds or crosslinks. All forces between the chains, except at the points of crosslinkage, are 
ignored and each chain is assumed to be free to take on any conformation in the confined region. 
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polymer solution, called a melt. 

For a system of M phantom chains, the partition function is written as a product of M 
Wiener path integrals (Section l2,lj) : 

For the sake of generality, each polymer i, defined by its path Rj(s) in (|4.1|) . has its associated 
contourlength Lj. For present purposes, this extra notation is unneccesary. Henceforth, all 
chains are chosen to be identical with respect to constitution and degree of polymerization, and 
will have equal length L. The notation C will refer to the total contourlength of all the chains, 
that is, C = Li- 

4.2 Confining the Melt 

The macromolecules are confined between two parallel planar surfaces or walls, as in Figure PTTl 
since symmetric problems are usually simpler to treat, and widen the scope of the problem. If 
one wall is situated at z = and the other at z = h z , then the problem of a network at a single 
surface is retrieved by taking the limit h z — » oo. 




z = h z 



z = 



Figure 4.1: Simple illustration of a very dense melt of phantom chains, confined between two 
walls. There is no tendency towards adsorption. 

The formalism chosen to represent the confinement, must ensure that all R (s) vanish beyond 
the boundaries. The partition function in (|4.1jl is modified, by noting that the path integral 
has to be carried out only in the allowable region [22]. In Section 12.21 a single chain was con- 
fined. Since the chains are independent, the partition function is the product of M single-chain 
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contributions 1)2.6(1 : 

M /•Ri(i)=R i 



Z {R,^;^} = AT TT / [PR* (*)] walls e-5^Jo HH ds (4.2) 

„ M 



x 



n 



e(Ri(s)-5)e(^-Ri(8)-i) 



(4.3) 



The normalization N refers to the number of configurations of a collection of M completely free 
polymer chains, each starting at R'j and ending at Rj in N steps, with Kuhn step-length I. The 
confining walls, which act as a constraint in the partition function, can be rewritten to resemble 
an attractive potential in terms of continuous arclength variables (Section 12. 2|) : 

A[Ri(s)} =ln[G(R(s)-z)]+ln[9(^-R(s)-z)], and 9 (R) = \ 1 lf ^ > (4.4) 

I otherwise 

Since the chains are only confined in the z-coordinate, the complete partition function for the 
confined melt will look as follows: 

M rR xl (L)=R x o r-L 

Z me i t = AAJ] / [VR xl (s)]^ V {-- R 2 xl {s)ds) 

"R yi {L) = Ry 3 fL . 

[VR yi {s)} exp - - / R 2 yi {s)ds\ 

RyM=R'y 1 M J® J 

Rzi(L)=R z , g rL 



X 



x / [VR zl (s)] exp - - / R 2 zi {s)ds 

Jr zz (0)=r' z M Jo 

M L 

+ ]T / A[R zi (s)]dsj . (4.5) 



It might be possible to approximate the infinitely deep well represented by In 0(R(s) ■ z) by a 
function that is more suitable to mathematical manipulations. 



4.3 Confined Network Formation 

Given the collection of phantom chains, they should be crosslinked to each other {polymer- 
polymer links) and to the confining surface {polymer-wall links) in order to constitute a polymer 
network linked to the walls, Figure l4~2l 

In practice, polymer-polymer 2 cross links are usually formed by joining two segments from 
different chains by means of sulphur vulcanisation, peroxide or radiation cross-linking j^I]. In 
this problem, polymer adsorption to a surface does not play a role whatsoever. Wall-links are 



2 Henceforth, polymer-polymer links will mostly be denoted by the term bulk-links, and polymer-wall links by 
wall-links. 
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Figure 4.2: A simplified illustration of a gel network formed between two parallel planar surfaces 
of width h z . The wall-links are situated an infinitesimal distance e from the wall surface. The 
dash- dot line is the line of symmetry at z = hh z . 

therefore synthesized in the same manner (for example irradiation) as bulk-links, but with a 
lower functionality. In a real cross-linking process, the links build up in time. One link will 
affect the neighbouring chain density, therefore influencing the position of the next link, and 
clusters of links may appear. As the crosslink density is increased there is an abrupt change of 
the melt from a viscous liquid to a solid, elastic gel that shows no tendency to flow. At this 
stage — the so-called gel-point — a giant cluster or network spans the whole sample. Theoretical 
explanations for how the transition actually takes place, still remain difficult |23| I24j . We shall 
only be concerned with the case of sufficiently high crosslink density; not a system close to the 
gel point. 

Here we assume that we are dealing with a homogeneous material where local clustering of 
crosslinks are negligible. This is a reasonable assumption to make if the network is formed from 
a dense melt system, or the crosslinker density is high enough. Furthermore, it has been shown 
that macroscopic quantities, for example the free energy, are insensitive to microscopic details 
such as inhomogeneities in crosslink-density |51j . 

4.3.1 The Origin of disorder 

Both types of linkages are permanent, and impose certain topological constraints between the 
polymers, given that the crosslink density is high enough. The randomness of the link formation 
is the origin of the quenched disorder in the system (Section 13.2(1 . What is the probability of this 
disorder? We can imagine that the chains — constituting a given, random link — were just 
touching prior to the permanent linking. The probability of the disorder is then given by the 
same weight as that of thermal equilibrium of the melt prior to network formation. This idea 
of an instantaneous linking mechanism, was first proposed by Deam and Edwards 9 . 
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4.3.2 Linking Formalism 

If a crosslink joins arclength point s* on chain i with point s* on chain j, then we have the 
constraint, illustrated in Figure H^T b) . 



R i ( a ?)=R i (^) 



(4.6) 



which is unbreakable and stays fixed at the same place on the chain. Mathematically, we use 
a Dirac-delta formalism to pick out the set of linkages, with the arc- locations of the crosslinks 
specifying the crosslink topology of the network. The crosslinks are constraints in the partition 
function of the system, which have to be approximated by a more tractable potential. 

We proceed, by considering a melt with chains touching at their future crosslink positions 
whilst incorporating the confinement formalism of Section f4.21 A simplified way to think about 
the manner of introducing crosslinkages, is by first adding iV c sliding links [2j in the bulk and 
N w touching links at both walls, at random. This is of course completely hypothetical, but will 
facilitate statistical calculations of fast crosslinking and accommodate the chosen probability of 
disorder in Section 14.3.11 





(a) 



(b) 



Figure 4.3: Two polymer chains, Rj(sj) and Hj(sj), in close proximity (a) modelled by a 
slip-link, and (b) a permanent crosslinkage constraint (|4.6|) . after a fast crosslinking procedure. 



After the touching and sliding links reach equilibrium with the system, they have a free 
energy F' given by 



. M 



1=1 

M L L 

x \yi / ds / ds' 5(Ri(s) -Kjis')) 
l i,j=i Jo Jo J 

x \ / ds j dxdy 5(Ri(s) — r)(x,y,0)) 



i=l 







M rL f Nv 
x| > / ds dxdy 5(Ri(s) - r](x,y,h z )) ,(4.7) 



1V1 

i=i 



o 



with T) a certain vector determining the positions of the crosslinks at the walls. The notation 
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[DRj(s)] ffi indicates that the path integration for the confined z-coordinates should only be 
carried out in the allowable region. Due to the z confinement the partition function (|4,7|) is not 
symmetric in space. Nevertheless, in each partition function, the crosslink arc-locations in the 
Dirac-delta functions are the same, since they are quenched variables. The crosslinks therefore 
act as a common combining factor between the confined and unconfined contributions to the 
partition function. 

4.3.3 Implementing Replicas 

The formation of the network on the walls is done by "freezing" the sliding and touching links 
(Figure f4.3|) of Section [4.3.21 such that the system reaches its final free energy F$, given by the 
logarithm of: 

e -F s ({x i} )/k B T = exp {-H a [{Ri}]/k B T) . (4.8) 
{Ri} 

In the above equation, the free energy F$ is that of a specific sample associated with a specific set 
of quenched variables S = {Xi}. The notation ^{r.;} denotes the integration over all possible 



conformations of the confined system with degrees of freedom frozen at S. In Section 13.2.21 it 
was shown that for a system with permanent constraints, the free energy of the system with 
a specific crosslink-state S should be evaluated first, followed by a disorder-average over all 
possible constraints. This procedure would give the correct experimental free energy T of the 
system: 

T = J V s ({X,}) F s [ {X^ } dX % = J V s ({Xt}) In Z [ ] dX t , (4.9) 

where Vs({Xi}) is the probability of arclength coordinates being frozen at a specific set {Xi}. 
Each sample's fabrication probability Vs — the initial probability of the chains when touching 
— can be be expressed as a formation Hamiltonian 3 . 

In Section T3, 2. 21 we averaged over the randomness of formation of a single bulk link by means 
of the replica method, replicating the system n times. This was done because the quenched 
average of a logarithm as in 1)4. 9j) is generally difficult to evaluate. We shall now extend this 
idea to a macroscopic number of permanent crosslinks. The probability Vs({Xi}) is now given 
by the ideas of Section f4. 3. II In the Deam and Edwards formulation the fabrication probability 
is manifested in the zeroth replica, coupled to the other n replicas (Hj. The free energy of the 
network system is then identified as the coefficient of n in the generalised partition function 
Z{n) of n + 1 replicas. Here, the generalised partition functions for the x and y coordinates 
coincide with the existing Deam and Edwards model for an unconfined phantom chain network 
[2]. In the z dimension, the walls act as a secondary localization of the chain coordinates. The 



3 The formation Hamiltonian is given by the a = part of the replica Hamiltonian 7i in 14.141 . 
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generalised partition function for the n + 1 replicas of a confined gel is given by 

» M ~ n M 2 

e -^w/%T =z(n) = Ary[ncR< o '(»)] w y[n ^II>R<")(»)] we -*I:,■i'( ft !■ , ) * 



1=1 



x i > i ds ds' n ^Rnsi-Rns')) 

i,j=l " u ^° a=0 

X [J2 L ds I dxd y II H^\s)-r,^(x,y,0)) 



i=i -° 



a=0 



x [ V / ds dxdy TT S(Kl a \s)-r,^(x,y,h z )) \ (4.10) 
L i=i 7 a=o J 

where ry^' determines the positions of the crosslinks at the walls. Each three-dimensional r/ Q ) 
is an element in a (n + l)-dimensional vector, given by 



rj = i rj(°) . . . rj( a ) . . . jy( n ) V 



(4.11) 



For the sake of easing the notation, we proceed by modeling the network by one large polymer 
of length C. 

In any network synthesized from many chains, there will be network "defects" |211 148 j . The 






(b) 



Figure 4.4: Illustrating three types of network imperfections: (a) A closed loop resulting from 
intramolecular cros slinking; (b) A dangling end (represented by the open circle) and (c) A trapped 
entanglement. 



first defect, in Figure l4~4T a). occurs due to the linkage of two points on the same chain, forming 
a closed loop that is not linked with any other chain. The second imperfection [Figure T4.4f b)] 
is when only one end of a chain is attached to the network, creating a so-called dangling end. 
Dangling ends and closed loops do not contribute to the elasticity or strength of a substance, and 
should be excluded when probing the elastic properties of the gel. In contrast to the previous 
defects, entanglements, Figure l4~lT c) . should not be omitted from calculations, since they can 
(in sufficient number) play the same role as crosslinkages (31 IE]. However, we shall assume the 
net contribution of above mentioned three defects to be negligible. 
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By rewriting the Dirac-delta constraints in Q4.1UJI as pole integrations 4 , we arrive at the most 
compact formulation: 



-Hn)/k B T = N £ N c \dii c f NJ.dfr 

2-Ki I 2iri 



n vr 



a=0 



x exp | - HI k B T - {N c + 1) log fi c - 2(N W + 1) log /x w | , 
with Tt giving the pseMcfo-Hamiltonian for network connectivity and confinement 

2 



(4.13) 



H/k B T 



^ ^ Jo 

a=0 JU 



C 



ds 



dR(°)(f 
ds 



+ A[RW(s)}+J2 ^[ r 



a=l 



+ / u c / ds ds' TT 5(R (a) (s) -R (a) (s')) 
Jo Jo a=0 

+ fi w ds dxdy Y[ ^(R (q) (s) -v (a) (x,y,0)) 

J° J a=0 

+ ^ J Q ds J dxdy \{ 5(^ a \s)-^ a \x,y,h z )) 



(4.14) 



where \i c and /i w are the chemical potentials of the bulk- and wall-links respectively. We have 
thus, in a convenient manner, represented the linking constraints as in a grand canonical ensem- 
ble. The crosslinks naturally confine the chain(s) to a well-defined region which must adhere to 
the defined region between the plates. 

The replicas 1 to n can describe strained versions of the O'th replica, and can therefore have 
different volumes and temperature. The integrations J in (|4.13|) imply integration over the n 
deformed systems, each with volume V = AV , with V the volume of the undeformed replica 
system and A the deformation tensor. After extension or compression of the system by A, 
the system will conserve its crosslinking topology that it had at fabrication, and therefore each 
replica must have the same set of crosslink arclength coordinates. The infinite deep square well 
potential representing the wall confinement (|4.4jl . is given by A in the unstrained system and 
by A for the strained replicas: 



i[R(s)] =ln[9(R(s)-z)]+ln[e(A 2 /i z -R(s)-z)], with Q(R] 



1 if R > 
otherwise 

(4.15) 



The constraint terms are exponentiated by implementing the following complex integral: 

B^ = ^_ I e^-^+V^dv, (4.12) 



where contour C encloses the origin. 
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4.4 Planning the variational calculation 

Only a handful of statistical mechanics models have been solved exactly. For the rest it is 
necessary to resort to some approximation method. The current replica model falls in the latter 
class due to the intractable path integrals in 1)4.13(1 . In this section we shall employ the Feynman 
variational principle, which follows from the inequality, 

(e x )>e< x > (4.16) 

which is valid for any real stochastic variable X, due to the concave nature of the e x graph |19j . 

4.4.1 Introducing transformation coordinates 

It is expected that the crosslinks will physically pin down the giant network chain in space, 
so that for any s, the coordinates of the different replicas, RA°)(s), RW(s), . . . R(") (s), will 
be correlated in some way. The coordinates of the undeformed (zeroth) replica, are free to be 
anywhere between the constraining walls. Before we show how to mimic a crosslink constraint, 
we introduce a new set of coordinates {X(°), X^ 1 ), v( m )} b m= i. ..(„_!)}, with X^ 1 ) a relative 
coordinate, and X^ ) being the centre-of-mass coordinate of all the replicas. The n— 1 remaining 
coordinates are simply rotations in replica space and give the deviation of the chains from the 
affine position. We define the new set of coordinates in the standard way |14j : 

x (o) = 5 +^=i A / fl i (4 . 17) 

1 (l + nAf)i 

(m) = 1 y e (2irima)/n R (a) ^ m = 1, 2, . . . (n - 1) (4.19) 

where the j's are Cartesian indices. These coordinates define an orthonormal transformation 
T with Jacobian equal to one. The A^-'s are the elements on the diagonal of the deformation 
tensor A. As a rule, a gel substance is only weakly susceptible to volume change during strain, as 
already mentioned in Section l2.4l Consequently, the A's are suitably defined by the isovolumetric 
deformation tensor A given by ()2.26() . The matrix form and other details of the transformation 
T are relegated to Appendix El 

4.4.2 Introducing a trial localization potential 

To make the problem more solvable in terms of the variational procedure, we introduce a trial 
potential to simulate the Dirac-delta crosslink constraints. Under strain the mean positions 
of the crosslinks may deform affinely, but it would be wrong to enforce the condition that all 
the crosslinks themselves should deform affinely. We therefore allow them freedom to oscillate 
around these affine positions. It seems reasonable to simulate the crosslink constraints by a trial 
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harmonic potential of the form 5 , first proposed in |14j : 

Q = h f * £ (4.20) 

1/0 i=x,y,z /3=1 

with % being the localization parameter and a measure of the limits within which each crosslink 
is allowed to fluctuate. If qi is small, the crosslinking is weak. The dimension of q^ 1 is length 2 . 
The inverse of qi defines the mean distance in which the crosslinks are localised in the i-direction. 
In other words, (qxQyQz)~^ defines the allowable volume which a crosslink may explore. In this 
chapter, we assume that the type of localization at the walls does not differ considerably from 
the bulk localization. Due to the symmetry of the problem in the x and y coordinates, we set 
Qx = Qyi and henceforth work with two scalar parameters q x and q z . 

5 Note that the replica indices in the transformed coordinates, {X' ' , X* 1 ' , Y^ m '} |{ m =i...(n-i)}i are now de- 
noted by j3, to avoid any misconceptions. Before the transformaton T, the formation (unstrained) replica was 
a = 0, and the remaining (a > 0) replicas were strained versions. Now, the replica (3 = represents the centre- 
of-gravity of all the replicas, with the remaining (/3 > 0) replicas being relative to it. Henceforth, any references 
to a indices should be understood as pertaining to the physical R a -° coordinate model. 
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4.4.3 Notation 



In this section we list all the constituents 6 of the model that we shall refer to frequently during 
the actual variational calculation. 



• The connectivity of the chains (approximated by one long polymer of length C) is denoted by 
the Wiener term W. 



n-1 



X< > 2 ( S )+X« 2 ( S )+^ \Y^(s)\' 



ds 



(4.21) 



The bulk crosslinks are denoted by /i c X c 

^ C X C = fi c [ C ds [ C ds' S(X^(s)-X^(s'))S(X^(s)~X^(s')) 
Jo Jo 

n-1 



I [ S(Y {m \s)-Y {m (s')) 



(4.22) 



m— 1 



• The wall crosslinks are denoted by /x w X w . [We assume the same chemical potential /x w for both 
walls. 1 



rC r+ca 

fi w X w = fj, w ds dxdy 

Jo J -co 



n-1 



5(xM(s)-vW( X ,y,e))5(xW(s)) f[ 5(Y^(s)) 



5(xW(s)-vM(x,y,h z -e)))8(xW(s)) R S(Y^(s)) 



(4.23) 



• The trial potential is 



6 



X 



n-1 

W 2 i \ " \v(™)\2 



711 — 1 



(4.24) 



• The wall-constraints are denoted by A. 

A = ( ds{A{T™xW +T™XV>) 







n n— 1 n 

A(TfX<P +T 1 /XV + T^WYM) I 
/3=1 m=l ' 



(4.25) 



• The constant chemical potential terms associated with the contour integrals, we group together 
as c, that is, 



c = -(JV C + 1) log/^ c - 2 (JV W + 1) log//. 



(4.26) 



6 The vector in l)4.23|l . defined in Appendix^] describes the wall-link locations, which were defined by the 
vector 7] prior to the transformation T. 
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In terms of the new variables {X^°\ Y^ m ^} and above notation 7 , the generalised partition 
function of 1)4.13)1 can be summarized as follows: 



e 



-F( Atw ,/ ic ,n)/fe B T = f ^ nt j g-W+^cXc+AtwXw+A+c = J j^j e -W/k B T_ (4.27) 



4.4.4 The variational Hamiltonian 

The motivation behind the variational (mean-field type) approximation is the relative ease with 
which one can obtain some understanding of the constrained system. 

The crosslinks play the role of a localization potential that restrains each chain to some mean 



path in the gel. If the confined gel is strained by A 1)2.26)1 the average trajectory of a Gaussian 
chain will also transform affinely. However, in the strained state (corresponding to replica 
indices a > 0) the chains will probably tend to deviate from their mean paths [3] , depending on 
the degree of crosslinking. The fluctuations should not be ignored, but rather simulated by a 
suitable trial potential. In Section [4.4.21 an harmonic trial potential 1)4.20)1 was elected, because 
it is simple enough so that the relevant path integrals in 1)4.13)1 can be calculated. 

We model the constrained system variationally by means of two variational parameters, 
namely the localization parameters q x and q z . Using the trial potential Q, the parameters can 
be found by minimizing the variational free energy J- VSir . Using the Feynman variational principle 
the expression 1)4.27)1 changes to 



e -F(At w ,/i c ,n)/fc B T > J ^ int j e (<Q+McX c +^ w X w }-W-Q+A 

[int] e (Q+McXc+MwX w } ( J g'^j 

= e ~ F LA<ixAz)/ k Ji T (4.28) 

where (. . .) means averaging with respect to the Green's function given by 

r. „ n—l 

J g' = J [OT(°»][DX«][ J] VY™] ^-w-Q+a . (4.29) 



771=1 



~' H var/ fc B T 



Calculating the Green's function in ()4.29)) seems intractable at this stage, especially due to the 
A term for the X^ 1 ) and Y( m ) variables. For now, we circumvent this problem by implementing 
another trial potential that only contains the centre-of-mass coordinate X( ), and which is 
defined as follows: 



c 

Aq = / ds In 



o 



G (r 00 x(°) (s)) Q(h z - T 00 xi°) (s)) } (4.30) 



7 The shortcut notation [int] refers to § § J [OX (0) ] [PX (1) ] [ U"~=i ^Y (m) ] 
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such that we now choose to have to work with an exactly solvable Green's function given by 

/n— 1 
[PXW][PXW][TT PYH] e-W-Q+Ap _ ( 431 ) 

"1=1 = e -Hvar/fe B T 

In terms of the partition function, it is required to calculate 

g] (4.32) 



Z(n,/i w ,/i c ) = ( e (McX c +MwX w +Q+A-Ao) )0 f 



> e <^cX c + Mw X w +Q+A-A )g Gj (4.33) 

= e -F va ,A<lx,qz)/k B T (4.34) 

The free energy of the network can now be performed by using the simple trial form (|4.31|) 
for the distribution of the chains. The resultant variational free energy is then minimised with 
respect to the trial function, to give the best possible approximation to the real free energy T 
(for the specific choice of Ti VSLT ): 

T « min { F var } . (4.35) 

The average in the exponent in (|4.33j) contains a reduced wall-constraint, A — Aq, which we show 
in Appendix El to be negligible, 

(A-Ao)~0, (4.36) 



under the conditions of our variational calculation and a softened wall potential A. 
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4.5 The variational calculation (for (A — Aq) = 0) 

The actual calculation of the variational free energy J r va , r (qx,Qz) involves two steps. The first is 
calculating the Green's functions needed in order to average the exponent in Q4.33JI . The second 
is evaluating the weighted averages. 

4.5.1 The Green's functions 

The task of computing the Green's function Q in Q4.31JI is equivalent to evaluating three different 
analytically solvable propagators 8 . 

For the x and y coordinates, the Green's function corresponding to the "centre-of-gravity" 
replica {[3 = 0), Q = f [int] e _w , satisfies the partial differential equation 

(Ts ~ I V 4 0) ) G ° = 6 " ^ (0)(SO) 6 {S ~ S>) ' (437) 

for a system in the transformed volume, namely a (n + l)-tuple cuboid, given by 

V = V JJ (l + nA?)* . (4.38) 

i=x,y,z 

The path integral for the remaining replicas (/3 > 0) includes the trial potential Q which 
mimics the fluctuations of the strained (a > 0) from the unstrained (a = 0) replicas. In this 
case, (|4.31j) simplifies to evaluating Q = f [int] e~ w ~^. The solution satisfies the differential 
equation 

(s-l^+f-^'w) M*V.xV:\>-J\)-'(xV-xV) *(»-»'). («») 

and is related to the propagator of a quantum particle 9 in an harmonic potential |19j . 

For the confined z coordinate, there is an added wall-constraint, which exclusively influ- 
ences 10 the centre-of-mass coordinate xi°^ . The Green's function, Q = J [int] e _w+A ° , for the 
"centre-of-gravity" {(3 = 0) replica satisfies a partial differential equation 

( §- g ~ | V^ 0) - A [X®(8)] ) SoiXtoX®; \s-s'\) = S - X®) 5(s - s>), (4.40) 
analogous to (|2.8jl . 

The different solutions to the various path integrals are well-known |53[ I19j , and are sum- 
marised as follows: 



Henceforth, the bottom indices shall refer to the integration variable. For example: Gm indicates to a path 
integration with respect to y/ m ^ , where m £ 1,2, . . . ,n — 1. The index i refers to the three Cartesian coordinates, 
except when otherwise stated. The original sample volume is V, and Vi the box length in the i coordinate. 

9 The kernel of a quantum particle of mass m in an harmonic oscillator of frequency lo can easily be retrieved 
from the polymer result 14.431 by the following substitution: s — > it, A — > S and q — > ^f-. 

10 This follows from our choice of a Green's function in (|4.3ipi . 
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^(I + tiA 2 )^ \27r£\s-s'\ 



1/2 



(x< >(.)-;r< V))' 



for i = {x,y}. (4.41) 



g (X| )( S ) ) X|°)( S '),| S - S , |)=^ | [MT°]exp{-Ajf X^ds} 

x exp { ^ In ie(T oa xW) - T 00 Af»)] ds} 



OO p_2-2 



Vl + n\ 2 h, 



E 



fa'?- 1 „ „'i 

TIT 1 1 



sm 





sin 




Vl + n\ 2 h z 







(4.42) 



27rsinh hlqi\s — s'| 



i _ Si [ X W 2 W + xf» 2 ( 5 ')]co s hI f<Il | s - s '|-2xW( )x( 1 '(«') 
g 2 sinh ^igjs-s'l 



(4.43) 



g m tf m) (s),Y< m \ S '), \s s'\) = TV / e-* l^"'l 2 «*»- Y./o £ iv/ m, l 2 - 

_3y( m ) (sjy*" 1 ) ( s ') 



27r sinh i^j|s — s'| 



1 „. [^ m ' 2 (,H^ m ' 2 ( ,' )|ln ,i,i(,. 

g 2 sinh 



(4.44) 



4.5.2 The Averages 

The average ( Q + /i c X c + /i w ^w ) in Q4.32|) in terms of the trial Hamiltonian 7i var , is evaluated by 
means of the Green's functions in Section The full calculation is presented in Appendix [O 
and here we shall simply list the results. 

The Bulk Crosslinks (/x c X c ) 



The weighted average of the polymer-polymer link contribution, defined in (|4,22|) and calculated 
in Appendix lC.il is given by 



(^ C X C 



fC rC 
He ' 



f dsf ds> J] S(X^(s)-X^(s'))5(X^(s)-X^\s')) 



n-l 



i=x,y,z 



m=l 



3/i c 



2 VI + nA 2 h z V 871-3 



2\ n/2 



£ 5 



3£ /£ 



A(l + n/A) 2vr^ ^ 



(4.45) 
(4.46) 



where ^4 is the cross-sectional area of the undeformed volume containing the sample. The chain 
cutoff -length £ c , over which the chain is not flexible but stiff, is of the order of magnitude of the 
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Kuhn length I. In order to compute the average (X cz )g of the zeroih replica, we limited the 
calculations to the limit of total chain length larger that the spacing between the plates, that is 

h z < Vcl of (trm 

The Wall Crosslinks (fi w X w ) 



The wall crosslinks are completely specified by the vectors i/(°\x, y, e) and u^°\x, y, h z — e), for 
crosslinks situated an infinitesimal distance e from the wall "surface" : 

v {0) (x,y,e) = (x VT + n/A, y y/l + n/X, e y/l + nA 2 ) (4.47) 
u^(x,y,h z -e) = (x y/1 + n/A, y y/T+ n/X, (h z - e) \A + ™A 2 ) . (4.48) 

The weighted average of the polymer- wall link contribution calculated in Appendix IC. 21 is 



2^ / ds / dxdy 6[X®{8) - 4 0) ) 5(X^(s) ~ 4 0) ) ^(^ 0) (^) " 



J-oo 

n— 1 

x^wj^n^i 1 ') n w n) )«( i i ra V(^ B) )) ( 4 - 49 ) 

n 



l«m w7 r 2 e 2 £ ^»/a f . M (|r)()) 



7T 



/if (l + nA 2 ) 1/2 
which is valid only if h z <C £ and -f* > 1. 
The harmonic trial potential (Q) 

The average over the harmonic trial potential consists of n identical terms fAppendix lC.3j) that 
differ only in terms of the variational parameters q x and q z . 



(Q)= Ej/ ik J, +Eif 12 

u=a;,i/,z m=l 



(4.51) 



= - 1,3+ -q z £coth-q z £+ —q x £coth-q x £J (4.52) 

= ^3 + ^ + 2^)) for ^>1. (4.53) 

4.5.3 The variational Free energy 

The variational free energy (|4.34|) of n + 1 replica systems is given by 

iWn, q x , qz) = e'^n^/k^T (454) 

= e <McX c +^ w X w +Q+A-A > g f f g) (4.55) 



x — — e 6(i+nA?)h| 11 v e ~ — (2fe+9z) 4>56 
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In the above expression the integrals over the chemical potentials p c and p w can be evaluated 
using (|4.26j) and the identity (|4.12|) such that the variational generalised partition function 
can be written as follows: 



(X C )^(X W ) 



[$- — (qz+2q x )\ ^ j e 6 (l+nA|)h| 9 (4-57) 




£ 2 



3£ 



^ IL \A + nA? 27r ^* >/T + nA| 



In 



7T Z 



n(JVc-l) 



8vr 3 



i ^ 2 2/" / 2 \ n /2 

16/x w 7rV£ f qzQx 



h\ (l + nA 2 ) 1/2 



7T" 



e 4 1 



Jg 6(l + nA|)h 



7T- 



(4.58) 



In Deam and Edwards [0] the replica limit n — > is taken prior to minimizing the free energy 
with respect to q. The variational free energy ^^(qx, q z ) of the original gel system can be 
identified as the coefficient of n as follows: 



dZ var /dn 



n=0 



>2var(n = 0) 



(4.59) 



1 



2 (1 + c/p) 



E 



cA 



i=x,y,z 



DA 



+ ^ + 2 fe) - ^ + ^-^ ln^ 



(iV c + iV w - 1) 



V2vr 



(4.60) 



where 



2^ ln f and P = C,V - 



(4.61) 



Next the free energy should be minimized with respect to q x and q z to find the best approximation 
(for the given trial Hamiltonian function in (|4.31j) ) to the real free energy of the system. The 
resultant stationary points are isotropic, and deformation-independent: 



6 (N w + N c ) 
££ 



(4.62) 



Substituting the ^-values in (|4.6()|) we find that the free energy on deformation, J 7 , has the 
following upper bound: 



T _ k B T 



1 



2 (1 + c/p) 



E + 



cA 



+ 



(A^ 



DA 



6M 



+ 



3N n 



i=x,y,z 

(N c + N w -l) 



E ln 



3(N W + iV c 
vr££ 



(4.63) 
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4.6 Results 



In this chapter we extended the Deam-Edwards theory [2] to a polymer network that has formed 
at and between two confining, parallel plane surfaces. 

is 



The free energy on deformation, including only the Aj-dependent terms in (j£ 



1 



2 (1 + c/p) 



E ^ + 



i=x,y,z 



P 



+ 



6 hi 



(4.64) 



where N c and N w are the total number of bulk-links and wall-links respectively. The above result 
is only valid when the crosslink density is sufficiently high, such that we are definitely dealing 
with a solid gel. Secondly, the chain-density of the melt, and the effective contour length C of 
the polymer network contributing to the elasticity should be greater than the spacing between 
the wall h z . 

In this model the fluctuations of the crosslinks were simulated with an harmonic potential 

and variational parameters q x and q z . The physical significance of these parameters is that 

_i 

they play the role of localizing the fluctuations a distance proportional to q x 2 in the x and 

_i 

y dimensions, and q z 2 in the z dimension. The fluctuation volume that each crosslink may 
explore can thus be estimated as 



-3/2 



6(iV c + iV w 



(4.65) 



since the localisation parameters ()4.62(l are isotropic and A- independent. 

The above result is similar to the Deam and Edwards result jHj for an unconfined phantom 
network, with no wall-links. In the case of a free, phantom network the localization parameter 
is given by the above result (|4.62|) with no wall-links (N w = 0), that is qi = T£ ■ The Deam 
and Edwards result for the free energy on deformation 



•Ttje < k B T 



with 



1 N r 



2 (1 + c/p) 



i=x,y,z 



2 1 



(—, ~ 



and p = C/V , 



(4.66) 



(4.67) 



is consistent with the confined-network result Q4.64JI in the limit of a gel only confined by a single 
wall (h z — * oo), and having no wall-links (iV w = 0). 

From (|4.64|) the stress 11 /, in the specific case of a uniaxial, isovolumetric deformation (|2.26|) . 
is given by 

£tt 2 C 



1 dT _ k B T 
! V d\ V 



_^_( A _J, + £ Ai + n „ 



A 



(4.68) 



xl The stress is the elastic force per unit area of the undeformed cross-section of the sample. The stress was also 
computed for simpler systems ( Sections 12.41 and 12 . 5|l . 
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The underlined part in (|4,68|) coincides with the well-known result of the classical theory of 
high elasticity of polymer networks jSJ, apart from the different front factor, which may be 
attributed to wasted loops [5]. 

The wall-confinement introduces new terms in the free energy of the network, dependent on 
the spacing h z between the walls. For example, the typical confinement term , is consistent 
with earlier results for a stitch network (Section 13, 4|) and a single chain (Section 12.4,1(1 with an 
effective contourlength L. 

To summarize: So far, we have applied the Deam and Edwards idea of homogeneous crosslink- 
ing and consequently obtained a similar deformation free energy (apart from confinement terms). 
It was found that the harmonic potential Q, controlling the fluctuations away from the mean 
(affine), depends on the inverse distance between (iV c + iV w ) crosslinks and not on the strain 
A. This means that the affine macroscopic deformation of the gel-plate system, results in a 
microscopic deformation of the chains, which consists of an affine contribution and a non- affine 
fluctuation contribution. 

Unfortunately the model employed in this chapter is too crude to differentiate between the 
types of crosslinking. This resulted in the wall-links and bulk-links being treated in the same 
manner. Intuitively the degree of localization, manifested in the values q x and q z , is expected 
to depend on the spatial coordinate. 

In the next chapter a network, formed from polymer brushes, is investigated. Within this 
framework, it is possible to distinguish, in a simple albeit concise manner, between the localiza- 
tion that each type of crosslink imposes. 



Chapter 5 

The Brush Network 



In this chapter we consider a network formed from two grafted polymer brushes. The term 
polymer brush was invented by de Gennes [Zj to describe an architecture in which polymer chains 
are terminally tethered to a surface at a high density This is the case when the separation k 
between anchors, is much smaller than the coil dimension, as shown in Figure IHTT1 



Figure 5.1: Simplified illustration of a brush, consisting of polymer chains attached to a flat, 
solid surface permanently, with average anchor distance k between chains. 

The anticipated behaviour (within the replica scheme) of a surface-attached, crosslinked 
polymer includes some kind of localization of the polymer. In Chapter 0] we employed the idea 
of Deam and Edwards [2], which treats this localization to be homogeneous and translation 
invariant, that is q x = q z = constant in the trial potential Q of Section 14.4.21 The average 
polymer network between planar surfaces is symmetric about z = \h z (Figure f4.2|l . but certainly 
not translationally invariant in the z-direction. Moreover, we expect the polymer to be localized 
differently near the surface than in the bulk region away from the surfaces 1 . This implies that 
the variational parameters qi should depend on the spatial height z(s), of a polymer segment 
s, between the plates. The only coordinate in the transformed replica system (Section 14.4. 
which represents physical position of the chain segments, is the "centre-of-mass" coordinate 
X(°). The remaining coordinates {X^jW} | 

{m=i...(n-i)} are simply relative coordinates. 
Mathematically, it is thus possible to distinguish between different localizations, by using a trial 



lr rhis is especially so if the distance between the plates becomes sufficiently large. There exists also the freedom 
of choosing the bulk and surface crosslink densities independently. 
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potential of the form 

/n 
dsq 2 {z(s))^2^ 2 , a>0, (5.1) 

with the localization parameter q a function of z(s) = z ■ X^ ) (s). Unfortunately, the path 
integration of an harmonic potential, which is dependent on the arclength coordinate s and a 
mixture of configurational coordinates of different replicas, seems intractable. 

However, in a brush network each stretched-out chain is only attached by one endpoint 
(s = 0) to the grafting surface. Consequently, we choose the average localization to be solely 
dependent on the arc distance s, and not the spatial position of the chain. It is with this 
simplification in mind, that we set out to investigate inhomogeneous localization, in terms of a 
brush architecture. 

5.1 Formation of the Brush Network 

In order to obtain a surface-attached polymer network, the chains of two polymer brushes should 
be sufficiently crosslinked. Each chain has the same molecular weight N (they are monodisperse) 
and contourlength L, such that L = Nl. Each brush consists of \N W phantom chains that are 
permanently attached via endpoints to a solid surface at very high density, that is k <C ViL, 
Figure 15.11 Attachment occurs in the absence of adsorption: the chains are not attracted to 
the wall, as in Figure !53T a). In practice, real chains will tend to stretch perpendicular to the 
surface due to overcrowding, until the loss in entropy balances the enthalpic gain of segment 
dilution. We again choose to work with phantom (Gaussian) chains, which are not subject to 
excluded volume effects and therefore overcrowding should not affect them. 

The best way to realize a brush network, is to grow chains directly at the surface of the wall. 
This approach is called grafting and is shown in Figure IS~21 




Figure 5.2: Schematic (cross-section) illustration of polymers grafting from a solid surface. 
The black circles represent the initiator molecules. 

Firstly the surface is covered with a monolayer of initiator molecules, from which long 
molecules grow like a lawn 2 . Secondly, the network is formed, by crosslinking the chains. There 
are a multitude of chemical fabrication techniques that can be used to effect crosslinking |54j . For 



2 We assume, as in previous chapters, that the polymer chains are linked an infinitesimal distance e away from 
the wall surface. In this case, the e can be interpreted as the size of the initiator molecules, from which the chains 
are grafted. 
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example, if the molecules have photoreceptive groups, the chains are crosslinked by irradiation 
with UV light |4()1 130 j . The formation of a brush network thus differs from the one considered in 
Chapter because surface attachment and crosslinking are not performed simultaneously, but 
in subsequent steps. 

5.1.1 Confinement formalism 

Prior to network formation, we have two brushes, that is, a dense system of N w independent 
chains, already attached to the wall surfaces. The planar surfaces are brought within a distance 
h z of each other, so as to confine the grafted polymers in the z direction perpendicular to the 
surface. Mathematically, the confinement is represented by an infinite square well potential A, 
defined in (|4.4|) and depicted in Figure l5~3T bV 



A ads [K(s)' 



(a) 



A[R(s); 



(b) 



Figure 5.3: Characteristic plots of the spatial profiles of two different situations: (a) A po- 
tential well corresponding to a hard, adsorbing wall; (b) An infinite square well potential, for 
confinement in the z direction between two confining plates of width h z . 

The partition function for a collection of iV w confined, grafted chains is given by: 

Z {{R h K' t ;L}) =M\Y\ [2?Ri(s)W4 e~^=^° H^H ds (5.2) 

I r=i M(o)=r^ J 

iV w r 

= AT jj / [©RiWle-Zo^a^W-^^W)]*, (5.3) 

where AA" 1 refers to the statistical weight of iV w free polymer chains, each starting at R^ (on 
the wall) and ending at Rj in L/£ steps. 

5.1.2 Crosslink formalism 

After confining the polymer brushes, the chains should be crosslinked to each other, in order 
to constitute a surface-attached network, Figure 15.41 Let the set of crosslink locations on the 
chains be denoted by {S}. These links are formed randomly, such that the crosslink positions 
{S a } of gel sample a will differ from {Sb} of sample b. However, once the links are formed, they 
are assumed to be permanent within a specific sample. 
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Figure 5.4: An illustrated, simplified view of a brush network. The walls are a distance, h z , 
apart and the long polymer chains are fixed at e from the surface. 



A cross linkage that joins points s* and on chains labeled i and j respectively, is described 
mathematically by a Dirac delta constraint, 5 [Rj(s*) — Rj(s|)] . A product of N c 5-constraints 
stipulates the topology of a network having iV c crosslinks 3 . Since the set {5} of crosslinkages 
is unknown prior to network formation, we assume that the linkage positions are distributed 
according to a probability distribution V(S). Let the distribution of the crosslinkages in the 
resultant network, be given by the probability of the crosslinkages an instant before linking, 
that is, 

Ni 



i=l 



X 



n^[Ri«)-Ri(4)i 



{5} 



Bulk 

N„/2 

X 

i=l 



U S[H i (s = 0)-rj(x i ,y i ,e)] S[Ki(s = 0) - rj^ylK - e)] , (5.4) 
' 1 - ' ' 



Wall at z=e Wall at z=h z 



where the 77's are wall-link position vectors [see Section f4.3.2| . This simplified history- dependent 
situation, first proposed by Edwards , is equivalent to a system where the future crosslinks 



3 A more suitable notation for s* mentioned above, would be sfj, where i and j are the chains involved in the 
linkage, and k the degeneracy of the pair. For example: A small network, consisting of three chains, could have 
the following incidence set {5} = {(sji, sfi), (513,821), (s?2,sli), (823,832), (si 3 ,S3i) }• 

*Js(«s) 

Ri(s? 2 ) = R 2 (s|i) w .O^; Ra(«2s) = 

R 2 (4)\~~j J 

Ri(si) 
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are just touching prior to an instantaneous linking 4 . 
5.1.3 Replica formalism 

Since the set {S} is fixed for a specific sample, the set and its linking history will remain unaltered 
if the specimen is deformed in some way. Let -F({5}) be the free energy of a particular sample 
that has been subject to strain A (Section I2.4|) . after network formation. The effective free 
energy of elasticity, .F(A), is obtained by taking the quenched average of the sample free energy 
over all possible realizations of the arclength locations {5}, 

F(A) = J V({S}) F({S}) dSi = J V({S}) In Z({S}) dS { . (5.5) 

The averaging of the logarithm in 1)5.5(1 is facilitated by using the mathematical identity (j.3.12() , 
Z n = 1 + nln Z + . . ., and replicating the system n times. At this point, it is possible to 
incorporate the probability distribution P({S}) of {S}, when the chains are touching before 
strain, as a zeroth replica. The effective free energy can then be identified as the coefficient of 
n in the following generalized partition function [Hj: 

e -F(n)/k B T = J 2 (°)({S}) [Z({S})] n e* V {F Q /k B T) = Z{n). (5.6) 

The factor 5 , e Fo ^ ksT in Q5.6JI . is the normalization factor of the probability P(S), and the 
statistical weight of a system where the arclength variables are not constrained to be a fixed 
set {S}. Note that the crosslink points {S} are common to all the replicas, and so we can easily 
average over them. Since we are interested in a resulting network with a sufficiently high crosslink 
density, we assume that each long macromolecule is involved in numerous crosslinkages. The 
constraint — averaged over all possible positions of the links — that picks out the N c crosslinks, 
is given by 

Crosslink constraint = 

The wall- links are "replicated" in a similar manner [see 1)5.8(1 ] . 

We now proceed as in Chapter |1J by calculating the generalized Gibbs formula of n + 1 

4 In Chapter 0] the choice of probability distribution was done more comprehensively, but also amounts to an 
instantaneous "freezing" of links. 

5 The free energy Fo coincides with F' of a slipping link system in l|4.7fl . apart from different naming of constants 
and wall-link integrations. 



Ww iVw 

EE 
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replicas, 
Z(n) 



N„ n 



[flpRf (*)] [H Jj2?Bt a) ( S )]e-S«E&/o i [^*?(»)-^(H*)]* 

i=l ^ i=l a=l 

Nw f-L n - K r 

x [ E E / / ds ' II - R 5 a) ( s/ )] 

L i=l 3=1 J ° J ° a=0 

jV w „ n i 

i=l a=0 

x[E / dxd V II S[K^(0)-^ a \x,y,h z )) 



(5.8) 



where the notation J refers to integration in an unstrained zeroth replica system, and J to 
integration in the remaining n deformed replicas. The bulk crosslink points have the freedom to 
form anywhere along the length of the chains. However, the wall linking always takes place at 
s = for each chain, but the precise location on the xy-plane surfaces is unknown. Therefore, 
the wall-link constraint in (|5.8ft is written without an s-integration, in contrast with the situation 
previously formulated in (|4.7|) . 



The crosslink constraints are exponentiated by introducing chemical potentials (|4.12|) /i c and 
fi w for the bulk and wall linkages respectively. This leads to the most compact notation for the 
generalized partition function, 



Z{n) 



a=0 



x exp<^ - H(n)/k B T - iV c log/i c - iV w log// v 



(5.9) 



where 7i is a pseudo-Hamiltonian for network connectivity and confinement similar to the previ- 
ous Hamiltonian encountered in (|4.14|) . The next task at hand is to evaluate the path integrals 
in (|5.9j) . Thereafter, the free energy of deformation J-(A) is recovered by taking the the replica 
limit (ETT21 . 
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5.2 The Variational Approach 

Evaluation of the difficult path integrals in (|5.9j) will be tackled by the usual Feynman variational 
procedure together with the replica coordinate transformation (|A.1|) of Chapter 0] 

5.2.1 Trial potential for the brush network 

The reason for investigating a brush network in the first place, is the prospect of an inhomoge- 
neous, but integrable (trial) potential, that can simulate the Dirac-delta constraints. With the 
brush network we capitalize on the fact that it is the endpoints, at arc location s = 0, that are 
permanently fixed to the walls. We do this by introducing two constant localization parameters. 
The first, qo, describes the localization within a certain region near the wall attachment, and q\ 
the bulk region, as illustrated by Figure IB31 



— qo <7i ► 

o • e-^ s 

s = t L 

Figure 5.5: The s-line depicts the domains of the variational parameters, qo, q\, and r in 

(ton . 

It is now possible to simulate the cross-link constraints by the following trial harmonic 
potential [compare with (|4.2()j) ]. 



Q = h E ( [ ds fl8 x( " )2 + [ ds £ x(/?)2 ) > e [1, 2, . . . , n] 



(5.10) 



and thus model the system variationally by means of three variational parameters, qo, q\ and r. 
In ([5.10)1 the chain conformations were written in terms of more convenient, transformed coor- 
dinates {X(°),XW,YH} \ {m 

=l...(rc-i)}i as previously defined by the transformation matrix T 

in (fOjl . 
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In terms of the new set of coordinates, the generalized partition function can be summarized 6 
as follows: 



\[ / [int]e- w < + ^ x « 
i=i ^ 

n j m e 



In the above and future succinct expressions, we employ the following notation: 



(5.11) 
(5.12) 



The connectivity of the chains, each of length L, is denoted by the Wiener term 



W ( = - 



2( 



>(™)/ s )|2 



ds 



(5.13) 



The bulk crosslinks are denoted by /i c X, 



£/ dsf ds'6(xf>\s)-X?\s'))6(xP(s)-X?\s')) 
7=1 Jo Jo 

n-1 

J] <5(Y^( S )-Yf(.s')) (5.14) 



• The wall cross- links are denoted by /z w X w . The variables v stand for the transformed vectors 
describing the wall-linking [see also Appendix lA"| . 



dx dy 



n-l 



5(xf(0)-^(x,y,e))S(X^(0)) f[ <5(Y«(0)) 



+ *(xf (0)-^(x,y,h z -e)))S(X^(0)) J[ J(Yf»(0)) 



(5.15) 



• The trial potential is 



ql J ds 



xP + ]T|Y i( m)f 



•a? 



x (i )2 + g |Y ( m)|2 

m— 1 



(5.16) 



The wall-constraints are denoted by 



n n— 1 >. 

a— 1 ' 



(5.17) 



6 The shortcut notation [int] refers to § § J [<5X<°>] [5X«] [ n^l\ <5Y< m >]. 

The notation c represents the constant chemical potential terms associated with the contour integrals, that is, 
c = -(iV c + 1) log^ c - (7V W + 1) log^ w . 
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If a second trial potential Ao (|4.3U|) is introduced 7 , analogous to Section ll~5l we obtain the 
following generalized partition function: 

> TJ e <McX c +/iwX w +Q+A-Ao> g ( g\ (5.19) 



i=l 

-F var (n,5o,iji,T)/fcBT 



(5.20) 



where the notation (■ ■ -)g indicates averaging with respect to the Green's function Q. We will 
again make the Ansatz (|4.36|) . (A — Ao) — 0. The variational parameters qo, q\ and t are 
recovered by minimizing the variational free energy, F ww . 

5.2.2 Averaging 

In order to obtain the variational free energy F va _ T in (|5.20f) , the average ( Q + /i c ^c + /i w X w ) is 
evaluated in terms of the Green function, 

g = j [JX^pX^in SY {m) ] e- w ~® +Ao „ (5.21) 

m=l = e -H var /k B T 

which can be factored, depending on the replica index and Cartesian coordinate, into well-known 
Green's functions listed in Section r4.5l Since we are working with a model having two parameters 
dependent on s, the averages should be performed carefully for each region, for example: 



(/(-) 



J dX s dX T dX L g(X s ,X ;s,q )f(s) G(X T ,X s ;T-s,q ) g(X L ,X T] L-r, qi ) . f 

/ dX T dX L g{X T ,X ;r, go) G(X L ,X T ;L~t, qi ) II U < S < r 

(5.22) 

J dX T dX s dX L g(X T ,X ;T,qo) G(X 3 ,X T ;s-T,q ) f(s) g(X L ,X 3 ;L-s, qi ) . f T 

J dX T dX L g(X r ,X ;T,q ) g(X L ,X T ;L-r, qi ) 



This section contains the final results of the averages, but the full calculation is relegated to 
Appendix IdI 



7 Note that Ao is not really a trial potential, since it has no variational parameters like Q has. However, since 
it only depends on the "centre-of-mass" coordinate X' ', it is integrable. 
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The Bulk Cross-links (/x c 



The weighted average of the crosslink contribution, defined in (|5.14j) and calculated in Appendix 
ID. H is given by the following lengthy expression: 



( 

] . i 27V w r 2 



3/Xc 
2 



3n 



+ 



VRJl + nX* 



+ 



/gi_\¥ J 2N w (L-t) 



+ 



TTlh z ^l + n\ 2 z 
£ c - y/t* + {L- r) 2 



L + r y/% + (L - ^ - L + r - 

2 V4 2 + (i-r) 2 -L + r J 2 L ^ + (L _ r) 2 + L _ T - 



(5.23) 



+ 



9o Qi 



7r(go + gi) 



2 



AN w t(L - t) 



V n 7l + nA 2 * ih * V 1 + nX l 



4 + VW+L 2 - VW+^ 



-V lt + (L-rr -Lln[ ^ + { ±Z^- L+ ^] -r 
V c V ' L s/£l + I?-L J 



In [■ 



V^ + (£-t) 2 -L + t- 



where the constant £ c is the chain cut-off length 8 . Equation (|5.23j) is valid when the polymer 
brush chains are long, and when the relation, V N£ S> h z holds. 

The Wall Cross-links (/x w X w ) 



The brush chains are attached an infinitesimal distance e from the planar surface. It is shown 
in Appendix lAl that any wall-link, after the transformation T (|A.1|> . is specified by the following 
two vectors, 



i/ (0) (x, y, e) = (x \]\ + n/X, y y/l + n/X, e y/T+ nA 2 ) 



(5.24) 



(x, y, h z -e)= (x x/l + n/A, y y/l + n/X, (h z - e) \A + nA 2 ) , (5.25) 



The continuous s representation of a macromolecule, is only acceptable for distances larger than £ c . At 
distances smaller than i c w £, the chain will appear stiff due to actual monomer bonds. 
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for the original bottom and top surfaces respectively. The weighted average of the wall-links, 
calculated in Appendix ID. 2\ is 

/ p+oo r 

(/x w X w ) = U w / dxdy ^xf^-z/^^^e^^xf^O)) [] ^(Yj m) (0)) 
\ • / -°° L m=1 



n-l 

+*(Xf )) (0)-i/f J )( a ;,»,fc,- e )))tf(xJ 1 J(0)) J] *(*f°(0)) 



m=l 



4^ w 7T 2 e 2 



(5.26) 
(5.27) 



The fact that the wall links in a brush network only occur at s = 0, is portrayed in the above 
(/-independent constant average, in contrast with the result for a "normal" surface-attached 
network l|I3U)) . 

The harmonic trial potential ( Q ) 

The average of the harmonic trial potential for inhomogeneous localization is given by the 
following expression: 

.2 f-r _ n-l n 2 r L n-l 

(Q) 



nt 



2 f'T 2 P Li \ 

-fl Mx ( ? )2 + Y,ri m) \ 2 ) + -fl ^(xS 1)2 +Ei Y l m) i 2 ))( 5 - 28 ) 

(^2[{L - t) q\ + t ql \ sinh cosh ^-(L - r) 
+ £<?0<?i j sinh ^ ([L - r]gi + rq ) + sinh ^ ([L - r]qi - rq ) 
x (qo cosh -|^t cosh -^{L — r) + cosh —^-{L — r) cosh -|^t] , (5.29) 

When the hyperbolic functions in (|5.29j) are written as exponentials, the average simplifies to 

n£ 

(Q) - T [(qo-qi)T + qi L] , (5.30) 

which is valid when the conditions, qo > and gi > e (l-t) ' are sa tisfiecl . 

In the limit of homogeneous localization, that is go = <Zi> the average (|5.28|) reduces to the 
familiar result given in (|4.53|) . 



5.3 The variational Free energy 

In order to obtain an approximate free energy of the original gel system (Figure l5.4|) . we first have 
to find values of qo, q\ and r that minimize the variational free energy -F^ar- At these stationary 
points, denoted by {c/q, gj, r*}, the variational free energy is expected be a good upper bound 
for the real free energy of a replicated system. In the last step the replica limit, which amounts 
to n — * 0, should be taken. However, for polymer network theories it is generally accepted that 
the replica limit and the minimization procedure may be safely interchanged Ej . We shall 
approach the problem in the manner of Deam and Edwards, by first taking the replica limit to 
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obtain the variational free energy T of a single system, followed by minimizing the resultant free 
energy. 
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The variational free energy 1)5.20(1 in 3 (n + 1) replica space is obtained by assembling all the 
averages of Section 15.2,21 

Z var (n) = e - F ^(nmm,r)/k B T ^ 31) 

iV w / r \ 

= Y\ e <McXc+MwXw+(Q+A_Ao ^ I G) (5.32) 

n Aw r N c \dfi c Mc(Xc ) I NJdy w » W (X W ) ^[(go-gi)r+ gl LI 
f „JVc+i e f „iv w +i e e 

/ 2 / \ — \ 



x - — e . 

A/" n 3 re 

= n(X c ) % (X w )^e-Tfeo-9i)T+?ii] e - 6(1+^1)^ i£. (5.34) 

In Section 15.2.21 some of the averages were simplified in the limit y/Nl 3> h z . For a dense 
network in a relatively narrow confinement, it is expected that the localization near the walls 

should not be appreciably different from that further away in the "bulk" region. We portray 
this expectation 9 by letting 

«D = (l + 7)9i> (5.35) 
With the above substitution, the variational partition function becomes 

5.3.1 The Replica limit 

The variational free energy of the original three dimensional system, say can be identified as 
the coefficient of n in ([5.6)1 . by employing the following procedure 9 , 

- f ^ T - 8 lr ( r:'or - (5 - 37) 



This substitution is also chosen to simplify the mimimization task at hand. 
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which results in the following variational free energy for the brush network: 
F(qi,7,r)/k B T = 

m ™(T^ \ ILX 2 N W 7T 2 3 ^ 4(1+7) ^ 3 (qi 
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(5.38) 



At this point, all that remains to be done, is to minimize the above expression with respect to 
the variational parameters {<?i,7,t}. 



5.3.2 The minimization problem 

The formidable task of finding the stationary points {(/f , 7* , r* } , where T has a global minimum, 
is equivalent to solving the following equations simultaneously, 



dqi 
dr 






on region r G (£ c , L — £ c ) 



(5.39) 
(5.40) 
(5.41) 



together with the condition that the Hessian H |25| of 
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(5.42) 



(91:7,t)=(^,7*,t*) 



and any second derivative, for example ^p"|(g£,7*,T*)j are both positive at the stationary point. 
Solving the problem analytically seems intractable due to the transcendental equation (|5.41|) . 
However, it is possible to find solutions, q\ and 7*, for (|5.39|) and l|5.40j) . The localization 
parameters are found to be 



Qi 



6JYc 



and q$ 



6A^ c (l+ 7 *) 



(5.43) 
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where 7* is still a function of r. The extremum points of T are calculated numerically. Next, the 
value t*, which minimizes the free energy is found by plotting the free energy as a function 



of r, and identifying r* as the value, on region (£ c , L 



where J-{t) has a global minimum. 



5.4 Analysis and Results 



In this section we probe the physics of the system, without beforehand obtaining closed form so- 
lutions for the variational parameters {<7i,t, 7} that satisfy the minimization equations (|5.39j) — 



3D. 



5.4.1 The homogeneous brush limit 

In a homogeneous brush network the crosslinks impose a uniform localization on the system. 
This situation is obtained, by letting r = or 7 = in (|5.43j) and (|5.38[) respectively. In this 
limit we obtain the usual Deam and Edwards 9 localization parameters 



6iVc 
£L 



(5.44) 



and the following free energy of deformation 
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(5.45) 

(5.46) 
(5.47) 



In the last step Q5.46JI we show only the important A-dependent terms. In the limit of uniform 
localization, we obtain a deformation free energy closely resembling the result ()4.64|) of the 
previous chapter 10 . The free energy ^uniform has the typical Deam and Edwards deformation 
term 5 -^cE an d front factor ~ (1 + c/p)~ 1 [see (|4.66|) ]. apart from a different c-value which 
now also depends on h z . The second term in (|5.46|) represents the typical confinement of a 
macromolecule between planar surfaces, and is common to the gel described in Chapter 21 as 
well as the stitch network of Chapter 



3 Here, the total number of crosslinks in the sample is given by N W N C 
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5.4.2 The single chain limit 

In Section lS.ll we mentioned that a polymer brush is (by definition) a very dense array of polymer 
chains, that is iV w 3> 1. However, if the crosslink density in the bulk is sufficiently high, we find 
a global minima for the free energy T , even if we only have a brush network synthesized from 
one chain, as in Figure EHJ a )- When N w = 1, we have a single chain network with one wall-link. 
The effect of the wall linkages relative to the bulk linkages, is then expected to be negligible, so 
that we are essentially dealing with a confined bulk network. A bulk network can be thought of 
as one giant chain, having no wall-links. 

Q Nc-1000, Nw=l 
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Figure 5.6: Illustrating the case of only one wall-link, N w = 1. (a) The plot of the variational 
free energy of the brush network T versus r exhibits a global minimum at t* ~ 305. (b) Plot of 
the localization parameters, qo and q% (|5.43j) . versus s for the brush network (solid line), in com- 
parison with a uniformly localized bulk network (dashed line), (|5.44j) . Constant plot parameters: 
L = 500, h z = 20, £ = £ c = 1, p = = 0.4, A = 1. 

From Figure l5~BT b) we can infer that the localization of a bulk network is expected to be the 
average of the inhomogeneous brush network localization contributions, qo and q\. 
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5.4.3 Localization of the polymer brush network 

Next, we examine the graphs in Figure 15,71 to see how the linking density N c influences the 
localization of the polymer. Each row in Figure 15*771 corresponds to a different value of N c . The 
first observation that can be made, is that the chains in the "surface" region (qo) are localized 
to a greater extent than the chains away from the walls (qi). This follows from the fact that 7* 
is positive in the middle column of Figure 15.71 and so we have that q\ < qo from Q5.43JI . The 
measure of the fluctuations of the crosslinks from their affine positions, given by the q- values, 
is greater for a chain that is lightly crosslinked (N c < N w ) as shown in Figure 15.71 than for 
a network where N c > N w . Since 7* < 1, the q- values Q5.43JI are primarily influenced by the 
number of crosslinks N c . 




(g) (h) (i) 

Figure 5.7: Investigating the influence of the crosslink density N c on the localization of the 
polymer. The leftmost column of graphs depicts the variational free energy as a function of r. 
Global minima were found at t* ~ 291 in (a), t* ~ 292 in (d) and t* ~ 295 in (g). These values 
are represented by vertical lines in the second column. The value 7* is given by the intersection 
of the vertical line with the j(r)-curve. The rightmost column illustrates that qo > q±. Also, 
the localization of the polymer decreases as the number of bulk-links decreases (|5.43|) ■ Plots were 
rendered with the following kept constant: L = 500, h z = 20, A = 1, £ = £ c = 1, p = 0.4. 
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5.4.4 Confinement and the free energy 

When a macro-molecule is placed in a geometric confinement, the total number of its possible 
chain conformations is reduced. This leads to a decrease in the entropy, or equivalently, to an 
increase of the free energy. 

However, the case of a confined, surface-attached network is more involved, since the crosslink 
constraints also influence the network chain conformations. We find that the entropy of the brush 
network increases as the distance between the planar surfaces gets narrower, with the volume 
of the sample kept constant. This phenomenon is illustrated in Figure l5~%T a) . where we plot the 
variational free energy J-(t) for two values of the wall separation h z . 




Figure 5.8: Illustrations of the variational free energy J 7 as a function of r, for (a) different 
values of h z at constant volume V, and (b) different values of polymer chain density p = N y 
for the the gel. Constant plot parameters were: L = 500, N c = N w = 1000, £ = t c = 1 and 
A = 1. 

At a sufficiently high cross link density and wall separation, the chains in a brush network 
may be thought of as highly stretched strings. When h z decreases, the tension in the chains will 
decrease, thereby allowing more possible conformations to the network chains, and an associated 
increase in the entropy. 

In this model the network is not confined in the xy-plane. This means that the average size 
of the chain is smaller than the x and y dimensions of the volume containing the gel sample. 
An increase in the box volume is equivalent to a decrease in the polymer chain density. As a 
result the free energy (entropy) of the network decreases (increases) , when the polymer density 
decreases, as shown in Figure l5~HT b'). The difference in the two graphs, lies in the fact that the 
chains are only confined and fixed in the z-dimension. 
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5.4.5 Elasticity 

It is possible to calculate the stress-strain relationship for the polymer brush network, even 
though a closed-form of the free energy cannot be obtained. This is so, because the localization 
parameters qo and q\ and the parameter r do not depend on the affine deformation of the 
system, that is, they are X-independent. By differentiating the variational free energy in (|5,38|) 
with respect to A, one obtains the following stress-strain relationship: 



/(A) 



V 



N r N m f. _. c. 1 tt 2 £LK 



(l + c/p) 1 ' p j 3hl 



A-VV + -A -— A 



(5.48) 



where the polymer density p and the factor c are defined in (|5,47|) , 

The factor A — 1/A 2 gives the form of the characteristic stress-strain plot so typical of the 
classical theory of rubber elasticity _48 . However, the deformation factor is altered by the 
constant front factor g = (l + c/p) , analogous to the result of the previous chapter (|4.68|) . and 
the Deam and Edwards [0] result for an unconfined network. The above result is consistent with 
the stress-strain relation obtained in (|4.68|) . for a surface-attached network where the wall- links 
and bulk-links are formed simultaneously. In the latter network, the wall-links have the freedom 
to form anywhere along the length of the polymer chains, in contrast with a brush network. 
This difference is portrayed in the iV w A term in (|4.68j) . which is absent from the brush network 
stress-strain equation. 

As an example, we list numerical estimates of the front factor constants for the Deam and 
Edwards model (unconfined network), and the confined network models of the previous and 
present chapter: 





Unconfined 


Confined 


Deam Edwards 


Real Confined Network 


Brush Network 


c 


0.17 


0.31 


0.14 


p 


0.4 


0.4 


0.4 


front factory 


0.70 


0.56 


0.74 



Table 5.1: The g- factors give a rough estimate of the percentage fraction of crosslinks that are 
elastically- active. The estimates in the example were computed for a network fabricated from 
N w = 1000 chains of chain length L = 500, chain cut-off length £ c = 1 = i, and confinement 
h z = 20, where lengths are in units of t. 



The front factor alters the classic jSJ elastic constant of N w N c k B T , since the number of 
elastically-able crosslinks are reduced from N W N C to g N W N C . This phenomena can be attributed 
to the fact that there was no concise way of excluding network defects, like closed loops, from 
the start. Note that the A term containing the front factor should not be translated into the 
elastic modulus of the system. In the new theories there are clearly other A 2 -dependent terms 
due to the confinement and wall-links, which will also play a role in the elasticity of the system. 

The second, /^-dependent term in (|5.48f) coincides with the stress of a single chain, of length 
N W L, confined between walls of width h z . Furthermore, the term -A, seems to be a distinctive 
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feature of a surface attached network, with the constant c depending only on the method of 
network formation. 

In this chapter we introduced inhomogeneous crosslinking in the simplest manner, by means 
of a confined surface-attached network fabricated from two polymer brushes. We saw that the 
polymer is localized to a greater degree in the "surface" region. The surface region is defined 
by s < r*, which turned out to be larger than expected. However, this should be ascribed 
to the fact that the calculations are only valid for the case when we consider relatively large 
macromolecules, \fT7t 3> h z . The stress-strain relationship (which contains the macroscopic 
information of the network) is not affected by the inhomogeneity that was introduced in the 
model. 



Chapter 6 

Concluding Remarks 



In this thesis we described three different confined, surface-attached network models: the stitch 
network, the real confined network and lastly the brush network. 

The stitch network in Chapter 3 captured the essence of confinement, and laid the ground- 
work (together with Chapter 2) — in terms of the relevant Green's function approach — for 
the more sophisticated models in subsequent chapters. Above all, the stitch network served as a 
preliminary exercise in treating disorder. The origin of the quenched disorder was the random, 
albeit permanent, crosslinking in the system. It was shown in Chapter 3 that one can handle the 
random crosslinking mathematically by applying the replica method. The controversial prob- 
lems sometimes associated with spin glass systems, are not applicable to our simple polymer 
networks, and it was sufficient to employ only the replica-symmetric approximation. The stitch 
network, can be described as a sum of smaller star-polymers (or two stitches), with four arms 
emanating from fixed positions on the parallel walls. Each bulk-link favoured the arrangement 
corresponding to minimum tension of the chains, namely at half of the chain contour lengths. 
It is therefore not surprising that the stitch network lacked the essential (classical) features of a 
typical network free energy. 

In Chapter 4, we implemented the ideas of confinement and disorder of the first two chapters, 
to develop a more acceptable model of network formation. This was done by adapting an existing 
model of network formation of Deam and Edwards jHj, to incorporate two parallel, confining 
surfaces (walls) and random wall-linking. The network was formed from pre-existing, confined 
macromolecules, that were crosslinked simultaneously to each other (polymer-polymer links) and 
to the confining surface (polymer- wall links). Calculations were performed under the following 
assumptions: 

1. The network was fabricated from a sufficiently dense melt, such that excluded volume 
effects were ignored. Consequently, in the phantom network, the chains could intersect 
each other and the possibility of entanglements was neglected. 

2. The crosslink density was sufficiently high, such that the gel was solid and not near the 
sol-gel transition point. 

3. The distance between the confining walls, h z was smaller than y/CA, the effective size of 
the macroscopic network chain. 
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4. The infinite square potential, which represents the hard wall potential, was softened in the 
framework of a variational calculation. 

Since we employed an harmonic potential identical to that of Deam and Edwards, with isotropic, 
homogeneous localization, we found the localization parameters to be equal, q x = q z , and anal- 
ogous to that of an unconfined network model. They were namely strain independent, and 
proportional to the mean square radius of gyration of a chain piece between two junction points. 
In terms of localization, the wall-links were treated as bulk-links, q = 6 ( Nc ^ Nvr ') , where C is the 
effective contourlength of the giant network polymer, £ the Kuhn steplength, and N c and N w is 
the total number of bulk-links and wall- links, respectively. The quantity q~ l l 2 gives a measure 
of the fluctuation of the network chains from the mean affine deformation path. In other words, 
it defines the relative diameter of a tube in which each chain is confined due to the surrounding 
crosslink constraints. 

In Chapter 5 we investigated a surface-attached network, fabricated via an instantaneous 
crosslinking of two pre-existing polymer brushes. This architecture was specifically chosen to 
facilitate an inhomogeneous localization scheme. This scheme modeled the effect of the crosslink 
and wall constraints by two constant localization parameters, which were arclength-dependent. 
It was found that each chain is localized to a greater degree near the surface at which it is 
attached, than far away from its grafting surface. The great advantage of this scheme, is in 
the fact that the wall-links were not treated as bulk-links, but played a role in the arc length 
distance over which each constant localization takes effect. 

Common to both the brush network and the real confined network, was the general form 
(apart from constants) of the stress-strain relationship /(A). Under the influence of an uniaxial, 
isovolumetric deformation, we obtain the following expression: 



with c a constant dependent on the chain length and confinement, and p the polymer density. 
The underlined part gives the characteristic curve of classic rubber elasticity models: a type 
of Hooke's law for small A's and a sharp rise in slope for larger strains. However, this slope is 
characterized by a different elastic modulus than encountered in the usual unconfined models. 
There are two new additions to the result of Deam and Edwards for an unconfined phantom 
model. The first is the c/p term, due to wall-links and the second, /^-dependent term, due to 
the confinement. The denominator (1 + c/p) gives a measure of the elastically active crosslinks 
in the system, and was also encountered in the Deam and Edwards model. 

In all of the calculations, we assumed a phantom model. Thus, an obvious improvement 
of the model would be to include trapped entanglements and excluded volume effects. This 
extension could for example be based on a non-Gaussian tube model or slip-link model ^H] with 
non-affine deformation. 

The variational calculation is extremely difficult when the localization potential depends on 
the relevant spatial position, say Z(s), between the parallel plates. For this reason, we chose 
a simple arclength-dependence. However, it might still be possible to find a potential that is 
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tractable, and more representative of the problem. 

Another future consideration is the study of the effect of fillers in the network. These volume- 
giving components can influence the network's response to deformation, by reducing the internal 
stress. It is known that a filler, like carbon black, leads to the reinforcement of network structure 
|29| . meaning increased stiffness, elastic modulus, tensile strength and abrasion resistance. The 
study of these improved properties, in the context of confinement and surface-attachments, may 
have a range of possible applications, and would be a worthwhile research endeavour. 



Appendix A 

Coordinate transformation 



The coordinates (s) = T^ a (s) define an orthogonal transformation with Jacobian equal 
to one 1 . In matrix form the transformation T may look as follows: 



( (1 + nA 2 )^ Xjil + nXj)-^ Xjil + nXj)-^ 

^JnXj (l + nX])-l 




\ 



1 ( 2irima ) /n 1 ( 2irima )/n 

s/n y/n 



\ 



(A.l) 



/ 



This is the most symmetric transformation method Qj . but there also exist other transformations, 
notably in [2], where all entries are chosen to be real. The factor (1 + nA 2 ) 1 / 2 ensures the 
orthonormality of the transformation. 

The wall-crosslink position vectors rf°^ are transformed in the same way as the polymer 
chain coordinates R/ a ) by the transformation T (|A.1|) as follows: 



,(0) 



(m) 



(0),y-n ■> (a) 
'lj jt" l^g=l A ] 'lj 

(1 + nA 2 )^ 
(l + nA 2 )5 

1 n 

1 \ " (2nima)/n (a) 



(a) 



m = 1,2, 



(n-1) 



(A.2) 

(A.3) 
(A.4) 



a=l 



The above calculation can be done explicitly for each Cartesian coordinate. For example: 

v2 



,(o) 



/(I) 



(m) 



X + Ea=l A * x 



(l + nA 2 )2 



\/l + nA 2 x 



A x x 



(1 + nXl) 

5 ( 2izim )/n 



II 



2-wim )/n 



m = 1, 2, . . . (n — 1) . 



(A.5) 

(A.6) 
(A.7) 



1 Cartesian coordinates x, y, z are represented by j. The index m runs from 1 to n — 1; .'. m = 1 <4> /3 = 2. 
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All except the zeroth replica variables vanish during the T transformation. 

After substitution of the coordinate transformation (jA.lj) . the pseudo-Hamiltonian of the 
generalised partition function can be written as follows: 



21 



Xi 0)2 ds 



Jo U£j 



2( 



A [T^Xf}+Tl a X^ s 



ds 



n—l 



+ e I a ' r °° x z°s + T z a x ^ + e 
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n-l 
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/ j z V Z S z s 



m=l 



+ I ds [dxdylfl s[T^(x(°\s)-uP( X ,y,e))+T^(xP( S )-vi 1 \x,y, 
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+ E T i m+1) " ( y i m) ( fi ) -^z m) (x,y,e)) ] + []^ 0a (xf)(s) - uf\x,y,h z - e)) 



m=l 
-.la 



s) - uP(x,y,h z - e)) 



o=0 



n-l 

+ E T i m+1)Q ( y i m) to - w i m) v> h * - ^ 



(A.8) 
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(A.9) 



Only the 5-dimension was shown, since this is the part of the free energy, which embodies all 
the constraints at this stage. In (|A,9|) we applied the fact that z/^ = for ^ 0, as was seen 
in lfA~8l) . 
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The average (A — Aq) in the 
variational principle 

We have to calculate the following: 1 

n 

e A = Q ( T f3o X (p)j Q ( hg _ T po X tf)j Yl e(Tf Q X^) @(\ z h z - Tf* XW) (B.l) 

a=l 

and e A ° = (T 2 00 Xj® )Q(h g - T 2 00 xj® ) = Q(xj® )®{ah z - xj® ) . (B.2) 
Now, A is a very hard potential that we choose to soften, by writing 

6 (Z) &(h - Z) ~ e"*" C 1 " 6 ^) @ ( h - z » , (B.3) 
where <f> is a big enough scalar, and Z stands for the transformed variables such that 

e < A> > e -^(n+l-(&(T z 00 X<«) 6( h z -T z 00 X,W))-ES=i < © P?" ^) ©( \.h,-T§ a ^ ) )>)/ B 4 ) 

e -(Ao> _^ e +0 2 (i-<e(TrxW)e(^-T°oxi°))>) ^ (B 5) 

where (. . .) denotes averaging with respect to the Green's functions listed in Section f4.5.11 

Next, we employ the fact that the integration limits, for the "centre-of-mass" coordinate of 
all the replicas X^°\ keep the argument of the theta-function positive and in the allowed region: 



t/T+^R^ (l+n\ 2 z )h z -R {0) 



(G {Xf ] + n^AxW) e ( y/1 + n\ 2 z h z - xj® - ni\xP)) = 1 (for A-term) (B.6) 
and < 8(lf) 9 (y/l + n\\ h z - X^)) = 1 (for A -term). (B.7) 



Substituting the above results for the j3 = and a = terms in (jB,4|) and (jB.5|) . we obtain 

g( A) „ -<P («-ES=i < (^" Q ^ Q) ) e CA»fc.-2!f- ) (B _ g) 
e"< Ao > ~ 1 . (B.9) 



1 Note that there is a summation convention implied with respect to the /3's, and a = (T^ ) 1 = V'l + ^^i- 
Also, since the transformation T is orthogonal, we immediately wrote (T^ 1 )"* 3 as its transponent T z a . 
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The n remaining terms in the sum in (|B.8|) can be written in terms of inverse Fourier integrals, 
where the Fourier transform of a stepfunction is 

3[6(Z)] = -L(-+7n^ . (B.10) 



\/27T \L0 

Since the coordinates Xz\ for (3 > 0, are simply fluctuations relative to the "centre-of-mass-of- 
all-the-replicas" coordinate X\ , we may expand the Fourier expressions in powers ofX^: 

(e (t^xW) e {\ z h z - T pa x^) ) (b.ii) 
= (i + i- J i^ {e -i^x^( s) A JM^^-iS-xFwA (B12) 



i 

i + - 
q 



du) terms 



+ o 



1 



,2m-l 



doj terms 



m = l,2,... (B.13) 



In the above average only the even terms [X^\s)] 2m survive the averaging over the Gaussian 
Green's functions [see Section r4.5.1| . and each of these terms contribute a ^ _ ( 2m_1 ) factor, where 
m is a positive integer. 

We thus have the following result: 

(A — Ao) ^ <p 2 — ( ^-independent terms ) + <p 2 ( n, 9 _-. ) - terms (B.14) 

q z V qz J 

In this thesis, we assume a high crosslink density. This implies that the localization parameter 
q is large, so that terms of the order q~ l are negligable relative to other terms, ~ q and ~ lng, 
which play the dominant roles in minimizing the free energy in ()4.6U|) . On the other hand, 
the scalar <j) is also large in order to decently approximate the hard potential A, and will tend 
to make the 1/q terms important. However, since these path integrals are not tractable, we 
approximate A by a soft enough potential so as to make the ^-influence insignificant. This 
assumption is thus inherent in the Ansatz of (|4.36|) . The centre-of-mass coordinate X^ is the 
only transformed coordinate that represents physical position of the polymer chains; the other 
coordinates are simply relative to it. In the variational calculation, the (3 > (or t/>-terms) are 
therefore only expected to play a relatively insignificant role in localizing the network. 



Appendix C 

The First Variational Calculation 



We set out here to evaluate the average ( Q + /x c X c + p w H w ) of Section l4,5.2l bv using the Green's 
functions listed in Section f4.5.11 



C.l The Bulk Cross-links (> C X C ) 

The (3 = term for the constrained ^-dimension is given by: 
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The solution to the integral in 1)0.21) follows from the following result: 

/ , v . „ . 2.^ • ^ sin(a - c)X sin(a - 26 - c)X 
I dX sin al sin 6X sin cX - 



(C.2) 
(C.3) 



4(a 



(a - 26 - c) 



sin(a + 26- c)X sin(a — 26 + c)X sin(a + c)X sin(a + 2b + c)X 



where a 
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a + 26 - c) 
6 



(a - 26 + c) 



4(a + c) 



(a + 26 + c) 



(C.4) 



7Tp 



and X = X 



(0) 



Each of the above terms 



Vl+nA 2 h z ' \/l+nA 2 h z ' \/l+nA 2 h z 

can only contribute to an integration of X = to X = VT+ n\ 2 h z if the quantities in round 



brackets go to zero. This immediately rules out the last two terms given in (jC.4|) since p,r, p > 1. 
The most straightforward combination to choose, is when a = 6 = c for each term, which implies 
that p = q = p, so that the three exponential functions reduce to one, with exponent ^i P h i L. 
Also, in the limit \f~L~l ^> h z , the exponential function (|G2|) is dominated by the first term in 
the sum, that is, p = q = p = 1. Consequently, for e <C h z , sin (7^) cos (7^) — jp, which leads 
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to the result in (|C.2jl . 

The (3 = terms for the x and y coordinates are given by: 

( 5{xf\s) - Xf\s')))g (C.5) 

_ / dx dx s ,dx s dx c g {x s ,,x ,s')g (x s ,x s ,, \ a - s'\) 5(x s - x s ,)g (x c ,x s , \c - s\) 
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- + 



2ir£\s -s'\ 



(C.6) 



where for example X s in the above is shortcut notation that stands for X^°\s). 

The (3 G {1,2, ... ,n} terms are identical, apart from different (/-values. In particular, for 
X^ the average is given by: 

( S(X^(s)-X^(s'))} Sl (C.7) 
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(C.8) 



Each Green's function in (|C7|) can be written as a sum of exponential functions multiplied by 
a product of eigenfunctions as follows: 
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where the last step follows from the assumption that the lowest eigenf unction, 

^0(X.) = (i)'e-^ Wa W, 



(C.9) 
(CIO) 

(C.ll) 



dominates the sum in (|C.9|) . This is the case when £qi\s — s'\ in Q4.43|l is large, so that 
cosh \lqi\s — s'\ — > ^expj^ls — s'\ \. In this limit, the integration in (|C.7|) can be repeated 
with the approximate Green's function £/i in ()C.10|) . such that the average becomes 



Since the Green's functions |s-s'|) of (1031) and g m (Y^> , Y^™ \s-s'\) of gSD 

for z E z}, all have the same form, the remaining averages are: 
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Putting all the averages together, the bulk crosslink contribution (|4,22|) becomes 
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C.2 The Wall Crosslinks (/i w X, 
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The wall cross-links are completely specified by the vectors u^°\x, y, e) and u^°'(x, y, h z — e), for 
cross-links situated an infinitesimal distance e from the wall surface. Letting a = \/T-j- nA 2 , the 
centre-of-mass average is given by: 
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The /3 = terms for the x and y coordinates are given by: 
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Originally, in terms of the R" coordinates, there was a clear cut distinction between the 
undeformed (a = 0) and deformed (a = 1, 2, ... n) replica s, and these were different by the 
deformation tensor A (|2.26|) . The wall constraint vectors were given by rj(x,y,h z — e) for the 
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a = replica, and 77 (-^, -M^, Xh z — e) for the a > deformed replica systems. Note that after 
the coordinate transformation T, the only vectors that do not vanish are 77' > for the top 

and bottom walls. When the network system undergoes an isovolumetric, affine deformation 
the averages are given by 
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Putting all the contributions together as dictated by (|4.49|) . results in the following: 
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The integral in ()C.31|) cannot be evaluated analytically. However, in the limit of £ > the 
integrand is constant, except at the boundaries near s = and s = £, as illustrated in Figure 
IC.ll This limit contains the localization factor q, which is found from (|4.62|) to be qi = 6 ^ N g£ Nc ^ 
for i = {x, y, z}. The magnitude of the non-constant contribution I err (near s = and s = £) is 
then of the order of X err ~ 2 (j\r +jv ) • This means that the assumption of a densely linked network 
(Xerr < 1) ensures that the integrand in ()C33|) can be safely approximated by a constant. The 
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Figure C.l: Plot illustrating the constant behaviour of the integrand in ()C.31I) and (|C.33|) . 

The solid line corresponds to ^£q > 1 (high crosslink density) and the dashed line to when the 
linking density decreases to ^£q <C 1. 



approximate solution to the integral 
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then leads to the result in (jC.32|) . It can also be obtained by rewriting the hyperbolic functions 
as exponentials and investigating the limit of large exponents. 



C.3 The harmonic trial potential 



Implementing only Green's functions Q\ and Q m which have the same form for all the coordinates, 
results in n identical terms. These solutions differ only in terms of the parameters q x and q z : 
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where i collectively denotes x and y coordinates. 
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The complete average (|4.24jl is found by integrating over the sum of the above terms: 
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Appendix D 

The Second Variational Calculation 



We set out here to evaluate the average ( Q + /x c X c + /x w X w ) of Section lB,2.2l bv using the Green's 
functions listed in Section f4.5.11 

D.l The Bulk Cross-links (/i c X c ) 

The (5 = term for the z-dimension corresponds with the average in the previous chapter (jC.3|) . 
and is given by 

{5{X<?H S )-X?\ S >))) go ^— =L=-, when p = r = p. (D.l) 

2v 1 + n\ z n z 



The solution to the integral in (|D.1|) follows from the result (|C.4|I and assumptions, previously presented 
inl(TTl 

The a = terms for the x and y coordinates are given by (where i and j are chain indices): 

( S(X^(s)-X^(s')))g (D.2) 
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J dX i0 dXiL Q (XiL , -^"io 7 L) 

i 
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(D.3) 



where V x denotes the length in the x-dimension of the original box volume, V = V x V y h z , containing the 
sample. The above result is for the special case when the polymer links with itself (i = j) to form a loop 
of length |s — s'\. When i ^ j, two different chains are crosslinked, and the average is s-independent: 

< S(xif(s)-X^(s')))g (D.4) 
_ J dXiodXjodXj S 'dXn J dXji J Qo(Xj S i;Xi( ) ,s)Qo(Xii J ,Xj S i; \L-s\)Qn(Xj S > , Xjq; s')Qo(XjL, Xj S i; \l-s'\ 
J dXiodXj dXiLdXjLQo(XiL, X i0 ; L)Qo(XjL, Xjo] L) 

(D.5) 



V x {l + nX 2 i y 

The above result occurs A^ w — 1 times. 

The rest of the bulk cross-link averages are performed in terms of Green's functions which are either 
qo or qi dependent, and therefore we have to employ the scheme in (|5.22|l . For large tq\s — s'\, we 
approximate the Green's function by Q\ given in IjC.lOf) . There are 3n remaining averages, given by 



79 



Chapter D — The Second Variational Calculation 



80 



(|D.6|) - (|D.9(I . where stands for any cartesian component of {X^ , Y| | m= i.. >T ,_i}: 
Qo Qi 



S(X^-X^)) =Jg (D.6) 

/ s<r, s' <r V Alt 



I s<t : s>>t y 7r (g + gi) 
^AT^-A^)) = JZ^KZ (D.8) 

= ,/f (D.9) 

/ s>r, s'>T V 27T 



Putting all the averages together, the complete bulk crosslink contribution (|5.23|) becomes 

.-. ( ^ C X C ) = ji c (X cx X cy X cz ) (D.10) 

21\ w h z 3 \ / go gi 
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V(l + n/X) 2%e\s-s'\J V2tt 

The above integrals should be evaluated carefully, since some of the terms diverge when s = 0. How- 
ever, the Gaussian chain model and continuous chain coordinates s, are only valid when we look at 
length scales larger than a certain length, say £ c . When s < £ c , the molecule is no longer a flexi- 
ble, continuous chain: it consists of monomers with rigid bonds. We therefore make the substitution: 
\s — — ► lim£ c ^o [ ( s — s') 2 + ]~ 1,/2 - The integrals are now straightforward to calculate and lead 
to the result shown in Q4.46J) . 



D.2 The Wall Crosslinks (/i w X w ) 

The wall cross-links are completely specified by the vectors v^(x, y, e) in l|5.24|l and v^(x, y, h z — e) in 
(|5.25(l . for cross- links situated an infinitesimal distance e from the wall surface. The wall crosslinks and 
endpoints (s = 0) of each chain will deform affinely together with the walls. Letting a = yl + nX 2 , the 
centre-of-mass average is given by: 

( S(X^ (Q) — z/°) (x, y, e)))g 

_ J dxi O) (0)dxi O) (£)g o (xi O) (0), Q e; 0) S(xj°\o) - ae) g Q (xj 0) (L), xj \o); L) 
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y/l + n\ 2 hl 



(D.12) 



since x| 0) (s = 0) = ae = y/l + n\ 2 . 

The x and y average for the "centre-of-mass-of-replicas" term, is 

< <5(Af } (a) i/<°> (x, y, e))) 0o = - ^ x2 , i = x, y , A, = Z= , (D.13) 

V^l + nXf vX 

where V is the volume of the original undeformed system. 
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The remaining averages are the identity because 

n-1 

(*(X«(0))) = <I[ «(Y W (0))) = 1, (D.14) 

m— 1 

because X (a) (s = 0) = Y (m) (s = 0) = (x, y,h z -e) = v {a) (x, y, e) = , for a > 
as shown in Appendix 1X1 

D.3 The harmonic trial potential (Q) 

Implementing only Green's functions Q\ and Q m which have the same form for all the coordinates, results 
in n identical terms for each of the s domains, for example for the coordinate: 

<Af )2 (*)) s<r (D.15) 
3g cosh £fs (q cosh £f(r - s) cosh \iq x {L - r) + qi sinh (t - s) sinh ^f{L - r)) 
q sinh cosh ^f-(L - r) + q x sinh ^(L - r) cosh l -fr 

{Xi 1)2 (s))s>r (D.16) 
3g cosh^f-(i- s) (q sinh l -fr sinh \i qi (s - t) + q x cosh ^(s-r) cosher) 
go sinh cosh ^f- (L - r) + 9l sinh l -f(L-r) cosh r 

The complete average ()5.16(l is found by integrating over the sum of the above terms, which is shown 
in the main text 
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